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' The massive field-theory approach for studying critical behavior in fixed space dimensions d < 4 

is extended to systems with surfaces. This enables one to study surface critical behavior directly 
in dimensions d < 4 without having to resort to the e expansion. The approach is elaborated for 

' the representative case of the semi-infinite l^l** n-vector model with a boundary term i cq jgy 0^ 

(~| , in the action. To make the theory uv finite in bulk dimensions 3 < d < 4, a renormalization of 

■ the surface enhancement co is required in addition to the standard mass renormalization. Adequate 
O ' normalization conditions for the renormalized theory are given. This theory involves two mass 
^ , parameter: the usual bulk 'mass' (inverse correlation length) m, and the renormalized surface 

I ' enhancement c. Thus the surface renormalization factors depend on the renormalized coupling 

I constant u and the ratio c/m. The special and ordinary surface transitions correspond to the limits 

m with c/m — > and c/m oo, respectively. It is shown that the surface-enhancement 
renormalization turns into an additive renormalization in the limit c/m — > oo. The renormalization 
, factors and exponent functions with c/m = and c/m = oo that are needed to determine the 

. surface critical exponents of the special and ordinary transitions are calculated to two-loop order. 

I ' The associated series expansions are analyzed by Pade-Borel summation techniques. The resulting 

, numerical estimates for the surface critical exponents are in good agreement with recent Monte Carlo 

^ ' simulations. This also holds for the surface crossover exponent "I>, for which we obtain the values 

O ' <l>(n = 0) ~ 0.52 and <I>(n=l) ~ 0.54 considerably lower than the previous e-expansion estimates. 

o : 

' \ PACS: 68.35.Rh, 05.70.Jk, ll.lO.Gh, 64.60.Fr 
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oo ' I- INTRODUCTION 

o : 

\ Sparked by the emergence of renormalization group (RG) methods at the beginning of the 1970s, the theory of bulk 
QQ ' critical phenomena has undergone a tremendous development in the past 25 years [0-^. Thanks to a very fruitful 

_ interaction with field theory, impressive progress has been achieved both in the theory of bulk critical behavior and 
"Y"^ . in field theory. While the latter has provided a rich variety of powerful tools such as Feynman-graph expansions and 

■ renormalized perturbation theory, on which analytical RG approaches could be based, the former has offered a wealth 
of challenging physical problems and served as a test laboratory for the application of new field-theory techniques. 

I ' One popular line of approach that has been extensively used with remarkable success are expansions about the upper 
' O \ critical dimension d* {— 4 for an ordinary bulk critical point) Q. The advantage of this technique is well known: 
The computational effort required for calculations to low orders in e = d* — d is relatively modest, in particular, if 
the simplifying features of such elegant schemes as dimensional regularization and minimal subtraction of poles |^ 
are fully exploited. As a consequence, the calculations can be — and have been [^|j7j — pushed to fairly high orders. 

A major reason for this computational simplicity is that the calculations can be performed directly for the critical 
(massless) theory. However, there is a price one must pay. The e expansion involves a double expansion in e and u, the 
^ renormalized coupling constant. In making this double expansion one by-passes the problem that the perturbation 
5^ ', series of the critical theory in terms of the massless propagator of the free theory is ill-defined for fixed d < 4 because 
of infrared singularities. In the dimensionally regularized theory, these singularities manifest themselves as poles at 
rational values of e which accumulate at d = d* as the order of perturbation theory increases |^,^ . Thus the problem 
of summing these infrared singularities arises. As stressed by Parisi pO|, without an additional hypothesis on the 
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summation of these singularities, any calculation based on the e expansion and the RG in this perturbative zero-mass 
scheme does not contain any information about the critical behavior in a fixed dimension d < d* . 

In practice, the e expansion often works amazingly well for critical exponents, even if truncated at order and 
extrapolated to d = 3 in the most naive fashion by setting e = 1. However, quantitatively accurate results require 
higher orders and sophisticated summation techniques P,|^,[rT|] . The extrapolation problem usually is more severe 
for other universal quantities such as amplitude ratios^ or scaling functions. One reason is that the results typically 
involve (e.g., geometric) factors or functions with an explicit dependence on d. Thus the question arises whether 
and which of these d-dependent terms should be expanded in e or rather be kept in the extrapolation procedure. 
As an empirical rule it has been advocated to choose the scale of u in such a fashion that a particular d-dependent 
geometrical factor is absorbed From a purely practical point of view, such recipes may well be useful. But they 
are hardly satisfactory since they neither have a firm theoretical basis we are aware of, nor ensure that all ambiguities 
of the extrapolation procedure are eliminated in a reliable fashion. 

The field-theoretic RG approach based on the e expansion has also been extended p^p^ to, and successfully 
used in, the study of critical behavior of systems with swr/acesQ JTsj-pOt . In the case of such systems an additional 
complication may arise: even at low orders of the loop expansion, the perturbative results may involve both geometric 
factors associated with the d dimensional bulk as well as others coming from the d—1 dimensional boundaries. Hence 
it may not even be clear how to apply the empirical rule just mentioned. 

From a fundamental point of view, approaches that work directly in a fixed dimension and therefore avoid the 
e expansion are clearly more attractive. An important one of this kind is the massive field-theory approach for 
fixed d < d* [ - 28 , ^ . Its merits are well known: Pushed to sufficiently high orders of perturbation theo ry a nd 
combined wit' 



1 sophisticated series summation techniques, it has produced values of bulk critical exponents 



with an accuracy comparable to that of the most precise ones obtained so far by alternative methods |p,r1JlllJ29,30[, as 



a, 23 1 



3^), we generalize the massive field-theory 



well as a set of amplitude ratios of barely inferior precision [g4 25 31|] . The method has also been utilized, albeit not 
to the same level of precision, to determine the universal ratio of correlation-length amplitudes for three-dimensional 
Ising systems in the analysis of critical behavior in various anisotropic and disordered systems p^-|35[|, partly 
even in general, non- integer dimensions 2 < d < 4 fs^ , as well as in studies of three-dimensional cjy^ theories describing 
the percolation transition and the Yang-Lee edge singularity problem j3j 

In the present paper (a brief account of which has been given in Ref. 
approach for fixed c? < d* to the study of critical behavior in semi-infinite systems. Such an extension is very 
desirable, both on account of the general conceptual reasons explained above, and for purely practical purposes. 
Recently extensive Monte Carlo calculations [^9|-|4^] have been performed for three-dimensional Ising models with 
free surfaces and for the adsorption of polymers on walls]^ ]45pqj39| , ^ . For most surface critical exponents these 
yielded values in reasonable agreement with the ones obtained by setting e = 1 in their e expansion to order [ p^ . 
For the surface crossover exponent $ however, the Monte Carlo estimates turned out to be 20-30% lower. 

These discrepancies were one of the motives for the present work. 

Our analysis is based on the semi-infinite n-vector model, which is the appropriate prototype model for studying 
surface effects on critical behavior [|l8|. In Sec. ||we briefly recall its definition and provide the necessary background. 
In Sec. IH we give normalization conditions for the massive field theory. Sections ^V|-[V^ are devoted to the analysis of 
the special transition. In Sec. |^ the general scheme of our approach is explained; then the Callan-Symanzik equations 
are given and utilized to derive the asymptotic scaling forms of the correlation functions near the multicritical point 
describing the special transition. After a brief discussion of some general features of perturbation theory, our two-loop 
results for the RG functions are presented in Sec. ^ These are utilized in Section VI to obtain numerical estimates 
for the values of the surface critical exponents of the special transition in three dimensions by means of Fade analyses 



and Pade-Borel summation techniques. The ordinary transition is treated in Sec. VII. Again, two-loop results are 



given and exploited to obtain Pade-Borel estimates of its surface critical exponents for d — 3. Concluding remarks 



are reserved for Sec. VIII. Various calculational details are relegated to Appendices A and B. 



II. BACKGROUND 



^For a review of results on universal amplitude ratios, see Ref. |12[ . 

■^For a review of surface critical phenomena, see Refs. and |20[ ; for an account of more recent results, see Ref. jl9| 
^ For background and references on polymer adsorption on walls, see Refs. and 
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A. The model 



Let = ((/'"(a;)) be an n-vector field defined on tlic half-space V = Rf = {x = {r, z) e R'^ \ r e R'^-^,z > 0} 
bounded by the plane z — 0, which we denote as dV. The semi-infinite n-vector model is defined by the Euclidean 
action 
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(2.1) 



Here TOq, uq, and cq are the bare mass, the bare coupling constant, and the bare surface enhancement^, respectively. 
Adding bulk and surface source terms to the action, we introduce the generating functional 



Z[J,Jj_;K, Ki] 



and the correlation functions 



Vcf) exp 



ay \ ^ 



(2.2) 
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(2.3) 



J=Ji=K=Ki=0 



For the functions G'^^'*^'°'°-' without 0^-insertions on or off the surface we use the notation G'-^'*^-'. The tensorial 
indices {aj^hk\ will be suppressed whenever no confusion is possible. The ultraviolet (uv) singularities of the theory 
should be assumed to be regularized by means of a large-momentum cutoff A. 

We shall also need the (bulk) analogs of these functions for the \(f)\'^ theory in the infinite space, i.e., with V — R'^. 
The easiest wa y to define these is the usual one where all boundary terms in the action (2.1) and the generating 
functional (|]^) are dropped, and periodic boundary conditions are chosen. We denote the so-defined bulk analog of 
G^^'^'^'^^dajj}; {X;}) as g\^{1^ {{xj}\{Xi\) and introduce their Fourier transforms g\^{1^ through 

GiT{{^jh{Xi)) 



(2.4) 



qi, ,Qj 



where the integral on the right-hand side indicates integrations J = J d'^{q/ 27r) over all c?-dimensional momenta 



, Qi . For the associated standard bulk vertex functions and their Fourier transforms we use the notation F 



(AT;/) 
bulk 



Qi, ■ ■ 

and f|,^ik\ respectively. 

In the case of our half-space geometry, where translational invariance is restricted to translations parallel to the sur- 
face, it is appropriate to perform Fourier transformations only with respect to (d— l)-dimensional parallel coordinates. 
We denote the {d — l)-dimensional parallel momenta associated with the operators <p{x ^ dV) and <ps = <P{x € dV) 
by lower case p's, and those associated with the insertions 0^ and 0^ by upper case P's. Parallel Fourier transforms 
are indicated by a hat; for example, the pair correlation function in this pz representation is written as G^^'^) (p; z, z'). 

Infinitely far away from the surface all properties must attain their bulk values. Hence the bulk functions gJ^^i']^'* 
can be obtained from G^'^'^'^'^'l by letting all N + I perpendicular coordinates Zj oo, keeping all relative coordinates 



Zk fixed: 



lim 



G(~."^^^°)(M;{z,}) 



(2.5) 



where {p} here stands for the set of all -|- / parallel momenta. 



^Upon mapping a semi-infinite (lattice) Ising model with ferromagnetic nearest-neighbor interactions of strength Ki between 
surface spins and of strength K elsewhere one finds that Co decreases as {Ki — K)/K increases For simplicity, we shall 

nevertheless use the term surface enhancement for co, rather than reserving it for (—Co) or (— co + const). 
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To proceed, it is necessary to recall a few well-known properties of the model ( J2.1D [ [18[ . Its phase diagram exhibits 
a disordered phase (SD/BD), a surface-ordered, bulk-disordered phase (SO/BD), and a surface-ordered, bulk-ordered 
phase (SO/BO), provided d exceeds the lower critical dimension dso/Buin) for the appearance of a SO/BD phase.^ 

The boundaries between these phases are the lines of surface, ordinary, and extraordinary transitions. They meet at 
a multicritical point, (mgjCo) — (mog,Cg^), called special point and representing the special transition. The ordinary 
and extraordinary transitions correspond respectively to the portions co > Cq^ and cq < c^J' of the line of bulk criticality 
TOq — ttiq^. The line of surface transitions separates the SD/BD from the SO/BD phase. At bulk criticality, we thus 
have three distinct transitions — the ordinary, special, and extraordinary transition. Of these only the ordinary and 
special one can be reached from the disordered phase. Since our present analysis is restricted to the disordered phase, 
only these latter two types of transitions will be considered. 

The restriction to the disordered phase simplifies the analysis considerably. One does not have to deal with a 
nonvanishing, and spatially varying, order-parameter profile {(f)[x)), and the free propagator in the pz representation 
takes the relatively simple form 



1 

2ko 



'K„\z-z'\ _ Co - ^0 g-Ko(z+2' 
Co + Ko 



(2.6) 



with 



60 . 



(2.7) 



The translation invariant first term on the right-hand s ide of ( |2.6| ) is the free bulk propagator. 

The perturbation series of the correlation functions (2^) in terms of the free propagator ( ^.6| ) can be regularized 
by setting G(p; z, z') — for \p\ > A. Whenever we do not use dimensional regularization, the theory is understood 
to be regularized in this fashion. 



B. Ultraviolet singularities for d < 4 



Let us first discuss the uv singularities of the theory. For bulk dimensions d = 4 — e < 4 the theory is super- 

renormalizable. Power counting shows that the uv singularities of the functions G*^^'^^) can be absorbed 
through a mass shift 

ml^ml + Srnl (2.8) 

and a surface-enhancement shift 

Co = Co -f Sco . (2.9) 

In order that the G^^^'^^^ be finite for 2 < d < 4 when expressed in terms of ttiq and cq, the contributions of order Uq 
to these shifts must behave as 

Sml - A^iuo/A')" (2.10) 

and 

Sco ~ A(uo/A^)'' (2.11) 

in the limit A — ^- oo. In contrast to Sitiq, which is known to be uv-divergent for d > 2, the shift Scq diverges only for 
d> 3. 



^This lower critical dimension is given by ciso/BD(l) = 2 and ciso/BD('^) = 3 for the Ising case, n — 1, and the n- vector case 
with n > 2 and 0(n) symmetry, respectively. In the presence of surface terms corresponding to an easy-axis spin anisotropy, 
a SO/BD phase is possible for n > 2 if d > 2. This case, studied elsewhere |^^, will not be considered here. The case d — 3 
with 0(2) symmetry of the Hamiltonian is special in that a surface phase with quasi-long-range order can appear, a problem 
which will also not be considered here. We shall also refrain from a discussion of {d = 2)-dimensional n-vector models with 
noninteger values of n in the range —2 < n < 2 ('loop models' see, e.g., El and its references). However, we shall estimate 
the surface critical exponents for the n case of polymer adsorption |45lM^, both for the ordinary and special transition in 
d = 3 dimensions. 
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C. Poles of the dimensionally regularized theory 



As is well-known , if the theory is regularized dimensionally, then the uv singularities of the bulk correlation 

functions manifest themselves as poles at e = 2/fc, with fc e N. These poles can be eliminated by means of an 
appropriate mass shift SmQ{e). We remind that the two- loop graph shown in Fig. 1 yields a contribution of the form 



ml-' hie), h{e)^ 



18 



1 



q,q 



, (g2 + 1) (^'2 + 1) (g + q')' + l 



(2.12) 



to 



f|,l(9 = 0) = X^;,k- (2-13) 
the inverse of the bulk susceptibihty. The integral has a simple pole at e = 1; i. e., /2(e) = ^2(e)/(e— 1), where Ri{e) 




FIG. 1. Two-loop Feynman diagram responsible for the poles in the dimensionally regularized 0* theory at d = 3 

is regular at e = 1 and whose value i?2(l) ~ l/327r^ can be calculated. Removal of the pole is achieved by 

n-f 2 1 



X 2^ \ 2/e 



18 327r2 (e - 1) 



Expressed in terms of ttiq and uq, the bare bulk functions rj^^j^j^ and G)^J^y are then finite at c? = 3. Yet, they also 
depend through logarithms on uq, and hence in a non-analytic fashion on it. 

Not only does this non-analytic behavior carry over to the correlation functions of our semi-infinite theory; owing 
to the appearance of additional uv (surface) singularities, it shows up already at one-loop order. To see this, consider 
the surface susceptibility 



(2.14) 



Xii(too,co) = G'^°'^^(p = 0; TOO, Co) 



(2.15) 



Its one-loop graph shown in Fig. 2 has a simple pole at e = 1, which can be removed by means of an appropriate 
choice of Jcq. 




FIG. 2. One-loop graph of the surface susceptibility xii having a simple pole at d = 3. The crossed circles denote points on 
the surface. 



We have 



n -|- 2 

Xii(too,co)"^ = Co + Too H — u^ra^^^ Si(e, cq/too) + e'(2-loops) 



(2.16) 



with 
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Ei(e,c)= Hdze-'^ f \ 

Jo Jp 2vV + l 

= Ji(e) + J2(e,c) + J3(e) , 

where = {1/271^-^ J d'^-^p, 

1 



C+ _^ 1 



Ji(6)e 



4Vr + l 
1 



2 Jp y/p^ + 1 1 + + 1 C + ^p^ + 1 



and 



^3(6) 



1 



p 4Vp2 + 1 1 + y'^^fl^ e - 1 



2-3^-i+./2 



r(e) 



r[(e + i)/2] 



(2.17) 

(2.18) 
(2.19) 

(2.20) 



The functions Ji and J2 are regular at e = 1; the above-mentioned pole results from J3. Upon expansion of the 
F-functions and computation of 



J2(1,C) = 



one arrives at the Laurent expansion 

Si(e,c) 

with 

Ri{c) 



cln2-ln(c+ 1) 
Stt (1 - c) 



+ Ri{c)+0{e-l) 



87r(e- 1) 
cln2-ln(c+l) 1 C^-lnTr 



87r(l-c) 

Demanding that the pole be canceled by Sco gives 



327r 



16n 



(5co(e) = -Mo 



i/e n + 2 1 



487r e - 1 



(2.21) 



(2.22) 



(2.23) 



(2.24) 



Expressed in terms of cq, mo, and uq, the bare susceptibility is finite at d = 3. At the level of our one-loop calculation, 
one finds 



n + 2 



■uo 



Ri[co/mo) - — In 

87r Uo 



0(2-loops) 



(2.25) 



The critical values 



ttiqc — SniQ and Cg = Cg — Scq of ttiq and cq pertaining to the special point would have 



to be determined from the conditions 



and 



Xii (^moc = \Jrnl^ + 5ml, 4"" = ^o'' + <^co j = 
The former is known to have the form 

rhlc - uj' M{€) . 

Similarly we have for the latter 

= uj" C{e) . 



(2.26) 
(2.27) 

(2.28) 
(2.29) 



The reason is that uq is the only dimensionful parameter remaining at the special point in the dimensionally regularized 
theory (with A = 00 and e > 0). In view of the non-analytic dependence of the susceptibilities on uq, it is clear that 
Symanzik's observation that Ai{e) cannot be determined perturbatively carries over to C(e). 

In the next section we describe an appropriate extension of the massive field theory approach that circumvents 
these difficulties. 
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III. NORMALIZATION CONDITIONS FOR THE MASSIVE FIELD THEORY 



Our aim is to study surface critical behavior at a bulk critical point. Therefore, a necessary property we ought to 
require from our approach is that the bulk critical behavior be treated appropriately. A convenient way to achieve 
this is to choose it in such a manner that it reduces for all bulk quantities to a well-established standard procedure. In 
our case this will be the conventional one based on normalization conditions (see, e. g., [p p^ , |2^ , |50| , |5]| ) . Alternatively, 
one could choose an approach based on minimal subtraction of poles of the massive theory, as described for the bulk 
case by Schloms and Dohm |29]. 



A. Bulk normalization conditions 

Starting from the bare bulk vertex functions r[j^[]^^(; m§, , uq), we perform a mass shift 

jtIq = + (5m^(e) 



(3.1) 



and introduce renormalization factors Z^{u), Z^2{u), Zu{u) (which are uv- finite for d < 4) as well as a renormalized 
dimensionless coupling constant u and renormalized fields such that 

(f) ^ [Z^{u)]^^'^ <p,cn , (f)^ ^ [Z^2{u)] UQ^Zu(u)m^u. (3.2) 

The mass shift and the renormalization factors are fixed through the standard normalization conditions 



bulk,] 



q=0 



(3.3a) 



j.g„(g; w, m) 



9=0 



(3.3b) 



q=Q=0 



(3.3c) 



fbulk.ren({9J;"'™) 



{q,=0} 



for the renormalized vertex functions 



Since the mass shift is sufficient to absorb the uv singularities of the bare functions r^^^j',^'' at d < 4, they become 
uv-finite when expressed in terms of m and u (or m and uq). The bulk renormalization factors can be written as 



(3.3d) 



(3.4) 



\^4>W)\ ^ = ^r[,^uik[9;"^o(™>'«)."o("i,w)] 



(3.5a) 



g=0 



\Z^i{u)Z^{u)\ ^ = f[,^^,fi^^[{0};mo(TO,'u),uo(m,u)] 



and 



[Zu{u)Zl{u)] ^ ^t^^l^^^[{0};ml{m,u),uo{m,u)]/uo{m,u) 



(3.5b) 



(3.5c) 
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B. Surface normalization conditions 



Consider now the cumulants G^^'*^) and G'^^^'^'^''^'^^ of our semi-infinite 0^ model. As we have seen, a sm^face- 
enhancement shift Sc is required to absorb uv singularities located on the surface. Hence we write 



Co = c + Sc 



(3.6) 



where c is a renormalized surface enhancement whose precise definition we still have to give. 

We also know that the surface operators 0s = 4>\z=o and ((/>s)^ should scale with scaling dimensions that are 
different from those of their bulk analogs 0(a;) and [(/)(a;)]^ with x ^ dV. This suggests the introduction of separate 
renormalization factors for these surface operators, which we do via the relations 



ll/2 



between the bare and renormalized operators. For the renormalized cumulants we thus have 
and 



(3.7) 



(3.8) 



(3.9) 



We wish to fix 6c and the new renormalization factors Zi and Z^2 by appropriate normalization conditions. To 
motivate our choice, let us recall the perturbation expansion of the momentum-dependent surface susceptibility 
Xii(p) = &^'^\p) to lowest order: 



G^°''^\p;mo,uo,co) = 



1 



Co + \/p^ + ™o 



-Oiuo) . 



(3.10) 



We choose the normalization conditions such that the associated renormalized susceptibility a nd its first derivatives 
with respect to agree at p = with the corresponding zero- loop expressions implied by ( [3.10 ), except for the 
replacements mo — > m and cq c. This gives 



G^°J\p;m,u,c) 



p=o m + c 



and 



dp'^ 



;m,u,c) 



p=o 2m(m -f c)^ 



(3.11a) 



(3.11b) 



The following condition fixes the normalization of insertions of the surface operator ■kcfy^, at zero external momentum: 



Gi",^''°''\p,P;m,u,c) 
This choice is motivated by the relation 

(5,(0,2;0,1)({0}) = 



p=P=0 



= (m + c) 



dcQ 



G("^2)(o) 



(3.11c) 



(3.12) 



Equation (3.11a) defines the required surface-enhancement shift Sc. Together with ( [3.3a ), it ensures that the special 
point is located at to = c = 0. The ordinary transition corresponds to the limit m — > at fixed c > 0. In this limit 
the renormalized surface susceptibility Xii,rcn ~^ . Hence the physical meaning of c is that of the inverse of Xii.ron 
at the transition]^ 



^Keeping c (and u) fixed while changing m requires, of course, that the bare quantities Co (and Uo) be varied with m. When 
exploiting the Callan-Symanzik equations below, we shall as usual hold these bare quantities fixed while varying m, so that the 
renormalized quantities u and c become m-dependent. 
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Equations ( 3.11b| ) and (3.11c) sp ecify the re norm alization factors Z\ and Z^i^ respectively, in a similar manner as 
the bulk normalization conditions ( 3.3b ) and ( 3.3c ) define Z^ and Z^i. The corresponding expressions are 



ZiZ^ ~ —2m(m, + c)^^^ G^*^'^^ [p; mo(m, u), uo(rn, u), co(c, m, u)] 



p=0 



(3.13a) 



and 



ZJ = -[ZiZ^] ^(m + c)^— G(°'^^ [0;TOo(m,u),uo(m,u),co] 

Cco 



Co— co(c,m,n) 



(3.13b) 



The above sets of normalization conditions ( [3.3a| )-( |3.3cj ) and ( 3.11a| )-(3.11c) define TOq, uq, Z^f,, Z^2, cq, Zi, and 



Z^2 as functions of m, u, c, and A. All Z factors have finite A ^ oo limits in the d < A case considered here. 

For simplicity, we consider the A = oo limit in the sequel. In our calculations described below we actually took 
A — oo from the outset, employing dimensional regularization. In this limit the bulk Z factors Z^, Z^2 , and Z„ become 
functions of the single dimensionless variable u. On the other hand, the above choice of normalization conditions 
( |3.11a )-(3.11c) implies that the surface Z factors Zi and Z^2 depend on both u and the dimensionless ratio c/m. 

In a full investigation of the crossover from the surface critical behavior characteristic of the special transition (for 
c/m <C 1) to that of the ordinary transition (for c/m ^ 1), it would be essential to carry along the dependence on 
the variable c/m. However, our main objective in the present work is the calculation of the surface critical exponents 
of the special and ordinary transitions. To this end, a study of the critical behavior in the asymptotic limits c/m — > 
and c/m — > oo is sufficient. As it turns out, there exist convenient procedures (see [ p8[ and below) which permit one 
to focus directly on these limits, avoiding the need to keep track of the detailed dependence on c/m. 



IV. SPECIAL TRANSITION 



Let us first consider the case of the special transition. In order to reach the corresponding multicritical point, we can 
safely set c = 0. This does not cause any problems in the theory provided the surface-enhancement renormalization 
has been performed. The desired critical behavior at the special transition can then be obtained by studying the 
massless limit of the resulting massive c = theory along lines analogous to those usually followed in the bulk case. 
It follows that the asymptotic critical behavior at this transition is described by the renormalized theory with the 
coupling constant u taken at m*, its value at the infrared-stable fixed point (and c set to zero). 



A. Normalization conditions at the multicritical point 



= the normalization conditions (3.11a)-(3.11c) simplify. The c = analog of (3.11a) fixes the critical value 
Expresse d in ter ms of renorm alized variables, it takes the form Cq = mf^{u) in the dimensionally regularized 



For c 
CgP of Co. 

theory. Equations (3.13a) and (3.13b) become 



ZlP{u)Z^{u) = -2m^-^ G(°'^ 



' [p; ma(m, u), Uo{m, u),cj'{m, u)] 



p=0 



(4.1a) 



id 

-m^ [Zl^{u)Z^{u)] — G''''^HO;mo(m, u),uo(m,u),co] 
Cco 



(4.1b) 



specifying renormalization factors Z^^(u) = Zi (u, c/m ~ 0) and Z'^ (u) — Z^2 (u, c/m — 0) appropriate for the analysis 

of the special transition. These renormalization factors enter the relations ( |3.8| ) between the bare and renormalized 
correlation functions G*^^'^^-* for c = 0, 



G(f„j^)(;m,.) = Z,-(^+^-^)/^(Z-)-^-^/^G(^^^^)(;mo,.o,c-) , 
and the corresponding c — analogs of the relations ( |3.9|) for G^^'^^'^'^^^. 



(4.2) 
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B. Callan-Symanzik equations 



By varying m at fixed uq and Cq , the Callan-Symanzik (CS) equations (of. Refs. m2^,^) of the correlation 



functions G^crl^Jp can be derived in a straightforward way. They read 

GL^-t?(;™,«) = AG.e„, 



with 



AG, 



where the integration is over the position X of the inserted cj)'^ operator. 
The RG functions appearing here are the usual bulk functions 



f3{u) — m 



_d_ 

dm 



and 



(u) = TO 



dm 



din Z^{u) 
du 



and the additional, surface-related function 



ri, in) = m— 



InZf (u) =/3(u) 



dlnZi"P(t 
du 



(4.3a) 



(4.3b) 



(4.4a) 



(4.4b) 



(4.4c) 



where |o indicates that the derivatives are taken at fixed bare coupling constant and surface enhancement (and cutoff 
A). ^ 

Just as in the bulk case, and as could be corroborated by means of a short-distance expansion, the right-hand side 
of (4.3a), AGrcn, should be negligible in the critical regime. The resulting homogeneous equations can be integrated 
in a standard fashion. 

In order to identify the crossover exponent $ we must also consider deviations Acq = cq — from the multicritical 
point. We use the expansion 



G(^'^'^)(;too,uo,co)=^ / .../ G(W/0(.^„,^„,4P) 

7=0 ^ JdvJdV 



(4.5) 



where the integrations /^^ are over the positions of the Ii inserted 4>1 operators. No infrared problems arise here 
because the massive theory is used. 

Expressing the right-hand side in terms of renormalized functions and the renormalized variable 



gives 



Hence 



Ac: 



Z:^{u) 



Acq 



Z;'^''''''\zf)-^'''-G^^^^''\-m,,u^^^^^ [ ...[ 

h-O V ^ 

Ii 



'-'rcn.sp ' (i ^) 



(4.6) 



(4.7) 



G(f„^)(; TO, u, Ac) ^ [Z,(u)]-(^+^)/^ [Zf (u)]-^/^ G(^'^)(; too, co) 



(4.8) 
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are well-defined renormalized functions.^ Since they depend on the additional dimensionless variable 

c = lS.c/m , 



their RG equations are analogous to, but differ from, (4.3a) through the replacement 



where 



m- \-f3{u) — 

am au 



rf^{u) = m 



dm 



d d d 

m— + I3{u)— - [1 + 77^P(w)] c — 
am ou oc 



(4.9) 



(4.10) 



(4.11) 



C. Scaling behavior near the multicritical point 



The CS equations given in the preceding subsection can be utilized in a familiar fashion to derive the asymptotic 
scaling forms of the correlation functions near the multicritical point. A detailed exposition of the derivation of scale 
invariance and universality of bulk vertex functions from the CS equations may be found, for example, in or 
elsewhere Since in the present case a completely analogous line of reasoning can be followed, we can be brief. 
In particular, we shall avoid carrying along the various non-universal constants (metric factors setting the scales 
of the relevant bulk and surface scaling fields), as would be necessary for an explicit derivation of four-scale-factor 
universality within the present massive RG framework. 

We shall need the familiar dependence 



(4.12) 



of the bare mass on r = (T — Tc)/Tc, valid for small deviations from its critical value m^^. We also recall that to, 
which is nothing else than the inverse of the (second-moment) bulk correlation length ^, behaves as 



near criticality, with 



where 77^2 denotes the value of the function 



m - (Too - TOoJ 



ri^2 (u) 



dr 



\tlZ^2{u) = (i{u) 



d In Z^2 (u) 
du 



at the infrared-st able ze ro u* of f3{u 



Integration of (4.4a), (4.4b), (4.4c), and (4.11) gives the asymptotic dependencies 

Z^-(u-w*)"/'"-to'', 



(4.13) 



(4.14) 



(4.15) 



(4.16) 



-(u-U*)''l''("'")/"-TO 



and 



{u-u*) 



(4.17) 



(4.18) 



^These should be distinguished from the previously defined c-dependent renormalized functions, which were related to the 
bare ones via c-dependent renormalization factors. 
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for u ^ u* or m ^ 0. As usual, u = P'{u*) and 77 = ri^{u*). Equation (4.18) can be substituted into (4.6) and (4.9) 
to obtain 



Ac ~ m-''°''("*' Acq , c - Acq . 

From the latter result we read off the scaling variable Acq with the crossover exponent 



(4.19) 



(4.20) 



Using these results, one easily sees that the CS equations of Sec. EVB yield the following asymptotic scaling forms 
of the correlation functions near the multicritical point: 



(4.21) 



Here /3 and (3^^ are standard bulk and surface exponents. The latter is related to the usual surface correlation exponent 
77|| via the scaling law 



where rj\\ in the present case of the special transition is given by 

sp . sp/ *\ 

■n{ ^■n + ViW ) ■ 



(4.22) 



(4.23) 



V. PERTURBATION THEORY 

We now turn to the explicit calculation of the surface renormalization factors Z^^ and and of their associated 
exponent functions rf^ and rf^ . In the one- loop approximation, this will be carried out for general space dimensions 
d < 4. However, in our two-loop calculations we shall restrict ourselves to the case d = 3. 



A. General Features 



The normalization conditions (3.11a), (3.111), and (3.11c) we have chosen to fix the surface counterterms all 
determine properties of the renormalized surface susceptibility Greif^ (p; m, u, c). For calculational purposes it is more 

(5.1a) 



convenient to express these conditions in terms of its bare inverse 1/G(°'2)(p). From ( |3lla| ) we find 

G(°'2)(0;mo,uo,co) 



= m + c 



Utilizing this, (3.11b) becomes 



_d_ 



p=0 



1 

2rn 



which shows that expression ( [3.13a| ) for Z\Z^ can be cast in the equivalent form 



(ZiZ^)-i = 2m^ \g^'^-^\v) 



Qp2 



m d ' 
= lim — ^ 

p=o p dp 



Likewise (3.13b) can be rewritten as 



(5.1b) 



(5.2a) 



(5.2b) 



It is useful to d ecompose the above inverse surface susceptibility into its free part, which is [G(p; 0, 0)] ^ = cq + kq 
according to (2.6), and a remainder due to perturbative corrections. Thus we write 
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(1) 



(2) 






(3) 



(4) 



FIG. 3. Feynman graphs to two-loop order of the no min ator g^Tg of the quantity d-o{p) in (5.5). Full lines denote the free 
propagator (p.q), dashed ones the reduced propagator (5.6). 



co + kq- ao{p) 



(5.3) 



To compute (Tq, we start from the following representation of the full propagator between two surface points in terms 
of E, the usual 'self-energy': 



\G 



\ S I s 



GTG\ 



T = E(1-GE)-^ 



(5.4) 



Here s| and |<; indicate that the left and right points arc located on the surface, respectively. A straightforward 
calculation yields 



CTo(p) = 



1+ s\G\sg'^Tg 

where ^ is a column vector whose components represent the reduced propagator 

g{p- z) = (co + Ko) G(p; z, 0) ^ e"""^ , 



(5.5) 



(5.6) 



and g^ is its transpose. 

The one-loop and two-loop contributions to ctq originating from the first two terms on the right-hand side of ( |5.5| ) 
are depicted in Fig. 3. Denoting the one from the graph labeled "(i)" by Ci{p), we have 



Co + no 



(5.7) 



where the term cx results from the last one in (5.5) 



B. One-Loop Approximation 

We now specialize to the case c ~ 0. Upon using the above results together with those presented in Appendix A, 
one can easily perform a one-loop calculation of the renormalization functions for general dimensions d < 4. Following 
Ref. , let us introduce a rescaled renormalized coupling constant ii through 

6 f47r)''/2 

u = bJd) u , bJd) = /, , ■ (5.8) 

^ ^ ' ^ ' n + 8 r(e/2) ^ ' 
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The advantage of this choice is that the expansion to order iP of the associated beta function /? = f3du/du takes the 
simple form 



/3(w) =-n{i-u)+o 



r,3 



Our results for the two surface renormalization factors of interest then read 

n + 2 u . ^^1 



Z^P = 1 



71 + 8 1 + e 



0[u 



and 



n + 2 1 



n + 8 1 + e 



. p /S-e 1 + e 3 + e.l 



u + 0(u') , 



where 2-F1 is the hypergeometric function [p4|. Substituting these power series into (4.4c) and (4.11) gives 



(5.9) 



(5.10a) 



(5.10b) 



(5.11a) 



n + 2 e 



n + 8 1 + e 



l-2^.^,f^,i±^;^;l 
' 2 2 2 2 



u + O 



(.») 



(5.11b) 



As a consistency check one can compute the pole parts (PP) of the Laurent expansion of the above renormalization 
factors at e = 0. One finds 



pp.^o [zT - 1] 



n + 2 u 
3e 167r2 



PP,=r 



Z 



sp 



(5.12a) 
(5.12b) 



As it should, this is in conformity with the one- loop terms (of the two- loop results) of Ref. ||T^, obtained by means of 
the usual scheme of minimal subtraction of poles at d = 4 for the massless theory in 4 — e dimensions.^ Accordingly 
we also recover the 0{u) expressions for the exponents functions 771 and rj,- of Ref. fl^ in the limit e ^ 0+: 



lim ri'l^{u) 

6-^0 + 



n + 2 u 



3 IGtt^ 



+ 0{u^) lim T]lP{u) + 0{u^) 



(5.13) 



From (5.9) one reads off the value u* = 1 in this one-loop approximation. Upon inserting this into (5.11), we get 

n + 2 e 



sp / * \ 

Vi [u ) 



n + 8 1 + e 



C'(2-loop) = ryjjP + C'(2-loop) 



and 



n + 2 e 



n+8 1 +e 



/3-. 1 + . 3 + e 1 
' " ' " ' 2 ' 2 



C'(2-loop) 



(5.14a) 



(5.14b) 



The reader may check that the Taylor expansion of these exponents to first order in e reproduces the known results 
(see Refs. [ p"6|jl^ a nd ref erence s therein). 
If we use ( ^.14a ) and ( 5.14b ) to estimate their values for c? = 3, we find 



77^P(n = 0) = 77^P(n = 0) = -i~-0.13 
77^P(n=l) = T]Y{n = l) = 



-0.17 



(5.15a) 
(5.15b) 



* Since in Refs. and a factor 2'*7r'^''^ was absorbed in the renormalized coupling constant, the quantity ti/lGvr^ here 
takes the place of the variable u of these references. 
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and 



^cP(" = 0)= g(l-41n2)~-0.22, 
r?^P(n=l) = i (1 -41n2) ~ -0.30 . 



(5.16a) 
(5.16b) 



Amazingly, the estimates (5.15) of this simple calculation turn out to be among the best ones for rj^ resulting from 
our much more involved two-loop calculations (see Tables | and On the other hand, our best two- loop estimates 
for ry^P differ appreciably from those listed in (5.16). 



C. Two-Loop Approximation 

At two- loop order we restrict ourselves to the case d ~ 3. Details of the calculation may be found in Appendix 
Here we just quote the final results. They read 



, , n + 2 , 12(n + 2) 
ZtZa, = 1 H u ■ -t;- 

^ * 2{n + 8) (n + 8)2 



n + 2 , , , , n + 14' 
A ^(1- In 2) In 2 



12 



48 



(5.17a) 



n + 2 f 1 

Z'%^1 + 21n2- - I u 

n + 8\ 2 



12 (n + 2) 
{n + 8)2 



„ n , n + 2, 9 2n+l 
A-B--ln2+^ln^2 + ^ 



u'^ + Oiu^) 



(5.17b) 



and 



n + 2 , 12 (n + 2) 



2(n- 



(n + i 



n + 2 , , „ n-10 

2A (1 - n 2 n 2 H 

6 ^ ' 48 



n + 2 
n + 8 



„ 1\ ^ 24 (n + 2) 

2 n 2 - - u - — 

2 n + 8 2 



A-B-^ln2 
2 



n + 2 , 2 17n + 22 
■ In^ 2 -I 



u'^+Oiu^) 



3 96 

where A and B are integrals originating from the two-loop graph (2) of Fig. 3 whose values 

A ~ 0.202428 

and 

B ~ 0.678061 

we have determined by numerical means (cf. Appendix ^). 



(5.18a) 



(5.18b) 

(5.19) 
(5.20) 



VI. SURFACE CRITICAL EXPONENTS OF THE SPECIAL TRANSITION 



We shall now discuss how the above perturbative results can be utilized to estimate the surface critical exponents 
of the speci al tra nsitio n. Our starting point are the series expansions of these expo nents i n powers of u, which are 
implied by (5.18a) and (5.18b). To generate these series, we substitute the expansion (5.18a) for rf^ into the following 
well-known scaling-law expressions for surface exponents: 



(6.1a) 
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(6.1b) 



/3i = 2 ^'^"^ + '"1) 



711 = ^ (1 - ^il) ' 
71 = 1^(2- ?7±) , 



^1 



A _ d + 2-T] 



(6.1c) 

(6. Id) 
(6.1e) 

(6.1f) 



^11 



cl- 



Ai 



/3i d - 2 + ; 



(6.1g) 



These scaling relations hold for the surface critical exponents of the ordinary transition as well; therefore, we have 
omitted the superscript "sp" . We also need the expansions of the bulk exponents v and rj. To the required order in 
it* , they read for the case d = 3: 



!/((i = 3, n) 



1 



n + 2 



2 4(n + 8) 



(n + 2)(27n-38) 
216(n + 8)2 



{u*f + ill*) 



and 



"^'^^^'")^ 27(n + 8)^ +'^[(-) 
We shall also consider the exponents 

af ^ a + ly - 1 + <P ^ 1 ~ iy[d- 2 - ??^p(u*)] 



and 



sp . 

a, 1 = a + 



2 + 2$== -i^[d-3-2r;^P(w*)] 



(6.2a) 



(6.2b) 



(6.3a) 



(6.3b) 



of the layer and local specific heats Ci(T) and Cii(T), respectively [Q. To obtain the expressions on the extreme 
right-hand side, we have substituted ( |4.20| ) for $ and used the hyperscaling relation a = 2 — dv. 
For each one of these surface exponents we arrive at an expansion of the form 



f{u*) 



k=0 



fk[u*t = h + hu* +f2{u*f 



O 



iu*f 



(6.4) 



As is known from the much studied bulk case (for background and references, see, e.g., Q), such series are asymptotic; 
they have zero radius of convergence. The reason for this is that the coefficients fk grow proportional to fc! as fc 00; 
more precisely, their large-fc behavior typically can be written as fk ~ Ck\k^~^A~^, where the factor fc! basically 
reflects the enormous growth of the number of diagrams contributing at a given order of the loop expansion. We expect 
that these features will carry over to the power series of surface quantities considered here. The large-order behavior 
of their coefficients and the values of the numbers A, b, and C should be obtainable by means of an appropriate 
extension of the instanton calculus utilized in the case of the \<p\'^ bulk theory. Furthermore, in view of the rigorously 
established Borel summability of the d = 3 dimensional model we may be confident that these series are 
Borel summable. 

In order to obtain meaningful numerical estimates from the above series expansions for surface critical exponents, 
we must invoke appropriate and sufficiently powerful summation techniques. The simplest procedure is to construct 
the table of Fade approximants [56| . This works well if successive elements Sn, Sn+i of the sequence of partial 
sums Sn{u*) = ^k=o fk{u*)'' vary little at low orders of N. A better and more sophisticated one is the Fade-Borel 
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method used in Ref. ||2l|. At the order of perturbation theory we are going to use it here, this involves the analytic 
continuation of the Borel transform 



OO ^ 

Elk (--^-.k 



/c=0 



(6.5) 



by a [1/1] Pade approximant. 

Our estimates given in Tables I-IV were produced as follows. For each exponent /, we rearranged the expansion 
as ///o = Mf = l + {h/h)u* + {hlfoWf or / + (1 - /o) = Mf = 1 + hu* + h {u*)\ depending on whether 
|/o| > 1 or |/o| < 1, respectively. Then Pade approximants of the type indicated in Tables I-IV were constructed 
for the so-defined modified quantities Mf, and [1/1] Pade approximants for their Borel transforms. For consistency 
reasons, these approximants were evaluated using the values of u*[d, n) one gets from the Pade-Borel resummed beta 
functions /3({2) at this two-loop order, namely]^ [pl| , ^ 



M*(d = 3,n = 0) = 1.632 



(6.6a) 



and 



u*{d = i,n = l) = 1.597. 



(6.6b) 



Finally, the resulting approximate values of the Mf were converted into estimates for the exponents by inverting the 
above equations defining Mf in terms of /. Note that we used the hyperscaling relations with d — 3. This is why 
our zeroth approximations (gathered in the column [0/0]) do not always reproduce the Landau (or e = 0) values 
0, 1, 0, i, 1, 3, 2 and 0, 0, — i, ^ of the exponents rf^ , . . . , i5^^ and 77*, . . . , $ listed in the first columns of Tables 
I/II and III/IV, respectively. 



TABLE I. Surface critical exponents of the special transition involving the RG function rff" for the case n = and d = 3. 





O1/O2 


Oii/02. 


[0/0] 


[1/0] 


[0/1] 


[2/0] 


[0/2] 


[1/1] 


R 




/(Ai.i^.r,) 


m 


-1.6 


-2.4 


0.00 


-0.204 


-0.169 


-0.079 


-0.107 


-0.126 


-0.134 


-0.130 


-0.133 


Ai 


-4.8 


-2.4 


0.75 


0.954 


1.006 


0.911 


0.886 


0.919 


0.921 


0.926 


0.921 


ri± 


-1.4 


-1.2 


0.00 


-0.102 


-0.114 


-0.027 


-0.017 


-0.059 


-0.063 


-0.064 


-0.053 


/3i 


0.0 


0.0 


0.25 


0.250 


0.250 


0.264 


0.264 








0.255 


711 


-3.6 


-2.1 


0.50 


0.704 


0.756 


0.647 


0.618 


0.660 


0.663 


0.668 


0.666 


71 


-4.5 


-2.1 


1.00 


1.255 


1.342 


1.199 


1.154 


1.209 


1.212 


1.220 


1.207 


61 


-2.3 


-1.6 


5.00 


6.020 


6.281 


5.576 


5.397 


5.711 


5.739 


5.764 


5.734 


Sii 


-2.4 


-1.5 


3.00 


3.816 


4.121 


3.481 


3.283 


3.578 


3.599 


3.628 


3.612 



TABLE 11. Surface critical exponents of the special transition involving the RG function rjl^ for the case n — 1 and d = 3. 





01/02 


Ou/02i 


[0/0] 


[1/0] 


[0/1] 


[2/0] 


[0/2] 


[1/1] 


R 




/(Ai,i^,»)) 




-2.1 


-4.8 


0.00 


-0.266 


-0.210 


-0.140 


-0.174 


-0.181 




-0.189 


-0.183 


-0.165 


Ai 


-12.7 


-2.9 


0.75 


1.016 


1.113 


0.995 


0.961 


0.997 




0.997 


1.006 


0.997 


V± 


-1.7 


-1.4 


0.00 


-0.133 


-0.154 


-0.056 


-0.040 


-0.084 




-0.089 


-0.091 


-0.067 


I3i 


0.0 


0.0 


0.25 


0.250 


0.250 


0.261 


0.262 










0.263 


711 


-8.2 


-2.6 


0.50 


0.766 


0.863 


0.734 


0.695 


0.737 




0.739 


0.747 


0.734 


71 


-11.0 


-2.4 


1.00 


1.333 


1.499 


1.302 


1.237 


1.305 




1.306 


1.320 


1.302 


Si 


-4.4 


-2.0 


5.00 


6.331 


6.814 


6.028 


5.779 


6.084 




6.101 


6.147 


5.951 


Sii 


-4.8 


-1.8 


3.00 


4.065 


4.650 


3.845 


3.553 


3.882 




3.894 


3.947 


3.791 



TABLE III. Surface critical exponents of the special transition involving the RG function rjl^ for the case n = and d - 





O1/O2 


Ou/02. 


[0/0] 


[1/0] 


[0/1] 


[2/0] 


[0/2] 


[1/1] 


R 




f{ai,v, ri) 


Vc 


-0.7 


-0.9 


0.00 


-0.362 


-0.266 


0.183 


0.055 


-0.144 




-0.168 


-0.160 


-0.119 


Ql 


-1.1 


-1.7 


0.50 


0.217 


0.280 


0.467 


0.399 


0.350 




0.336 


0.342 


0.342 


ail 


-0.8 


-1.1 


0.00 


-0.362 


-0.266 


0.109 


-0.021 


-0.157 




-0.180 


-0.170 


-0.140 


* 


-0.4 


-0.4 


0.50 


0.421 


0.427 


0.642 


0.657 


0.479 




0.474 


0.474 


0.518 



The n = value is given by the negative of the value of the fixed point denoted U in Ref. |57 
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TABLE IV. Surface critical exponents of the special transition involving the RG function r]l^ for the case n = 1 and d = 3. 





O1/O2 


Ou/02^ 


[0/0] 


[1/0] 


[0/1] 


[2/0] 


[0/2] 


[1/1] 


R 




f{ai,u, ri) 


n* 


-0.7 


-1.1 


0.00 


-0.472 


-0.321 


0.168 


-0.052 


-0.200 


-0.230 


-0.215 


-0.144 


ai 


-1.4 


-2.9 


0.50 


0.131 


0.230 


0.393 


0.304 


0.284 


0.268 


0.279 


0.279 


ail 


-0.9 


-1.6 


0.00 


-0.472 


-0.321 


0.042 


-0.153 


-0.226 


-0.253 


-0.237 


-0.182 




-0.4 


-0.4 


0.50 


0.397 


0.407 


0.649 


0.661 


0.470 


0.463 


0.463 


0.539 



The quantities Oi, O2, On, and 02i appearing in Tables I-IV are defined through the expansions Mf — 1 + Oi + 
O2 + . . . = [1 + On + 02i + . . .]^^ of the modified quantities Mf = ///o oi Mf = / + 1 — /o. Using the latter quantity 
to generate such truncated expansions in the case |/o| < 1 rather than simply factoring out /o yields better behaved 
'inverse expansions', i.e., series for MJ^. 

The bigger the absolute values of the ratios O1/O2 and Ou/02i are, the better is the quality of the resulting series 
for Mf and its reciprocal 1/Mf, respectively. All ratios Oi(i)/02(i) have negative sign or else vanish. Thus all series 
produced by these expansions are alternating, and hence adapted to the above-mentioned Pade-Borel summation 
technique. (If a series were not alternating, it would be unsuitable for this method because the [ 1/1 ] approximant of 
its Borel transform would have a pole on the positive real axis, i. e., inside the integration range |22|.) The estimates 
obtained via Pade-Borel resummation of the power series for Mf and l/Mf are listed in Tables I-IV as R and R~^, 
respectively. 

In most cases the resulting power series in u have second-order corrections 02{i) whose absolute values are smaller 
than those of their first-order ones. Thus the sequences of associated partial sums appear to be slowly convergent, to 
the available low order. Exceptions are some series involving ry^f, whose behavior is rather bad. In the first group of 
exponents, related to rjTi^) shown in Tables I and II, the most reliable estimates are obtained from the direct 
series for the exponent Ai, which appear to exhibit the best convergence properties. These estimates are Ai = 0.921 
for n = and Ai = 0.997 for n = 1. Substituting these along with the standard bulk values ||2^,^ v(n = 0) = 0.588, 
77(n = 0) = 0.027, i'{n = l) = 0.630, and 77(71 = 1) = 0.031 into the scaling laws (|6Ta|) -(|6li), we have computed the 
remaining seven exponents of this group. The resulting values /(Ai, v,?]) are presented in the last columns of Tables 
H and ||. By and large, the agreement with the results obtained from the individual expansions is quite reasonable. 
The deviations of the values /(Ai, v, rj) from the other resummation estimates might serve as a rough measure of the 
numerical accuracy achieved here. 

The situation is less favorable for the second group of exponents, rj* = 77^p(u*), whose estimates are given in 



Tables III and [V. Their series exhibit poor convergence properties. One should be cautious in relying on estimates 
derived from individual series expansions of an apparently divergent nature, like in the case of the crossover exponent 
<i>. In this group the exponent ai has the estimates with the least scattering. The best series convergence has Ui^, 
and the corresponding Pade-Borel estimates are ai — 0.342 for tt, = and ai = 0.279 for n = 1. Accepting these 
together with the bulk values of v and i] given above, the estimates of 7;*, an, and $ listed as /(ai, i', 7;) in the last 
columns of Tables III and |^ were derived via scaling laws. 

In order to check the quality of our procedure for estimating the numerical values of critical exponents, we have 
applied the same kind of analysis to the series expansions to second order of the bulk critical exponents. The results 
are shown in Table ^ and compared with the estimates of Ref. (last column) . The agreement gives us confidence 
in our method. 

The numerical values of surface critical exponents gathered in Tables ^-|V generally are in reasonable agreement 



both with previous estimates based on the e expansion as well as with those obtained by other means. For comparisons 
we refer to Section IIL8 of Ref. [Q, where e expansion estimates and estimates that had been gained by alternative 
techniques till 1985 are given, and to Table VI for more recent results. Note, however, that our estimates for the 
crossover exponent $ are definitely lower than the values <I>(7^ = 1) ~ 0.68 and — ~ 0.67 quoted in Ref. [ p^ . 
The latter were obtained by setting e = 1 in the e expansion of $ to order e^. On the other hand, recent IMonte 
Carlo simulations yielded the significantly lower estimates $(1) 0.461 ± 0.015 |||, $(0) = 0.530 ± 0.007 H, and 
$(0) = 0.496 ± 0.005 ||3|. Our present results $(1) ~ 0.54 and $(0) ~ 0.52 are fairly close to these values. 

To see whether comparatively small estimates for $(71, c? = 3) can be obtained from the e expansion, we have applied 
the analogous summation techniques to the series 

2$ = 1 + ai (71)6 + 02(71) , (6.7) 



whose coefficients are known to be |16 



77 + 2 1 i=0.125 for7. = 

2(n + 8) ] 1 ^0.167 for 71 = 1 
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TABLE V. Numerical estimates of bulk critical exponents. 





O1/O2 


Ou/02. 


[0/0] 


[1/0] 


[0/1] 


[2/0] 


[0/2] 


[1/1] 


R 




Rof. y 


n ^ 


u 


-7.0 


-4.1 


0.50 


0.602 


0.614 


0.587 


0.583 


0.589 


0.590 


0.591 


0.588 


1 


-4.9 


-2.5 


1.00 


1.204 


1.256 


1.162 


1.137 


1.169 


1.172 


1.177 


1.1615 


a 


-7.0 


6.2 


0.50 


0.194 


0.266 


0.238 


0.238 


0.232 


0.230 




0.236 





-44.1 


-13.6 


0.25 


0.301 


0.304 


0.300 


0.300 


0.300 


0.300 


0.300 


0.302 


A 


-6.0 


-2.7 


1.25 


1.505 


1.570 


1.462 


1.434 


1.468 


1.470 


1.477 


1.462 



n = 1 



V 


-27.7 


-5.9 


0.50 


0.633 


0.654 


0.628 


0.624 


0.628 


0.629 


0.630 


0.630 


7 


-11.3 


-2.8 


1.00 


1.266 


1.363 


1.243 


1.207 


1.244 


1.245 


1.254 


1.241 


Q 


-27.7 


2.8 


0.50 


0.101 


0.215 


0.115 


0.148 


0.115 


0.114 




0.110 


/3 


14.5 


405.8 


0.25 


0.317 


0.321 


0.321 


0.321 


0.321 






0.325 


A 


-17.5 


-3.1 


1.25 


1.583 


1.703 


1.564 


1.525 


1.565 


1.565 


1.575 


1.565 



n = 2 



V 


21.7 


-9.1 


0.50 


0.656 


0.685 


0.663 


0.661 


0.663 




0.664 


0.669 


7 


00 


-3.2 


1.00 


1.312 


1.453 


1.312 


1.273 


1.312 




1.322 


1.316 


Q 


21.7 


1.9 


0.50 


0.033 


0.181 


0.011 


0.086 


0.010 






-0.007 




7.2 


16.5 


0.25 


0.328 


0.334 


0.339 


0.340 


0.340 






0.3455 


A 


36.1 


-3.5 


1.25 


1.640 


1.816 


1.650 


1.609 


1.651 




1.662 


1.661 


n — 3 


V 


9.1 


-15.9 


0.50 


0.673 


0.709 


0.692 


0.693 


0.694 




0.695 


0.705 


7 


14.5 


-3.6 


1.00 


1.346 


1.528 


1.370 


1.333 


1.371 




1.382 


1.386 


Of 


9.1 


1.6 


0.50 


-0.019 


0.159 


-0.076 


0.042 


-0.083 






-0.115 


/3 


5.2 


9.5 


0.25 


0.336 


0.345 


0.353 


0.356 


0.357 






0.3645 


A 


10.7 


-4.0 


1.25 


1.682 


1.910 


1.723 


1.686 


1.727 




1.739 


1.751 



TABLE VL IVIonte Carlo estimates for surface critical exponents of the special transition in d = 3 dimensions. 



n = 



711 


0.805(15) 


Meirovitch & Livne, 1988 M 




0.714(6) 


Hegger & Grassberger, 1994 |43|] 


71 


1.304(6) 


Meirovitch & Livne, 1988 [^91 




1.230(2) 


Hegger & Grassberger, 1994 ]43|] 




0.530(7) 


Meirovitch & Livne, 1988 HZ 




0.496(4) 


Hegger & Grassberger, 1994143] 


n = 1 


/3i 


0.18(2) 


Landau & Binder, 1990 [ 


lOl 




0.22 


Vendruscolo et ai, 1992 






0.237(5) 


Ruge et ai. 1993 || 






0.2375(15) 


Ruge & Wagner 1995 ^4 




7ii 


0.96(9) 


Landau & Binder, 1990 [ 


40) 




0.788(1) 


Ruge & Wagner 1995 ^ 




71 


1.41(14) 


Landau & Binder, 1990 [ 


40) 




1.328(1) 


Ruge & Wagner 1995 |4 






0.59(4) 


Landau & Binder, 1990 [ 


40| 




0.74 


Vendruscolo et ai, 1992 






0.461(15) 


Ruge et ai, 1993 Q 
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and 



02 (n) 



87r2(n + 8) - (n^ + 35n + 156)] 



4(n + 8)3 

^TT^ - ~ 0.4645 for n = 0. 

16 256 

^7r2 - if ~ 0.5336 for n = 1. 



(6.9) 



The results are shown in Table VI. It is reassuring that the estimates obtained via Pade-Borel summation compare 
reasonably well both with our above ones based on the perturbation series at fixed d = 3 as well as with the Monte Carlo 
results mentioned. That these estimates deviate considerably from the values obtained from the [2/0] approximant 
( |6.7D at e = 1 seems to be due to the unusual largeness of the O(e^) term of In summary, we conclude that the 
values of the crossover exponent $(n,d) with n = 0, 1 and d = 3 are indeed significantly smaller than previously 
thought, being close to 0.5. 

TABLE VII. Estimates for ^{n, d — 3) based on the e expansion. 



n 


O1/O2 


Ou/02^ 


[0/0] 


[1/0] 


[0/1] 


[2/0] 


[0/2] 


[1/1] 


R 







-0.27 


-0.28 


0.5 


0.438 


0.444 


0.670 


0.740 


0.487 


0.482 


0.483 


1 


-0.31 


-0.33 


0.5 


0.417 


0.429 


0.683 


0.757 


0.480 


0.474 


0.475 



An interesting aspect of the above results is worth mentioning: We may be quite confident that the inequality 

al\ < (6.10) 

is satisfied for d = 3 and n — 1. For one thing, our best numeric al es timate based on the massive RG approach at 
fixed d = 3) is otiiiji = 1, 3) ~ —0.18. Second, the scaling relation ( |6.3b| ) can be rewritten at ci = 3 as 

af-^{n,3){n,d^3) ^ -2{v - ^) . (6.11) 

In view of th e various estimates given above it seems rather unlikely that $(1,3) will be larger than i/(l,3) ~ 0.630, 



so that ( 3.10 ) should be valid at d = 3 and n = 1, a conclusion whic h ma y also be reached for n = 0. 



This has important consequences. As has been shown in Ref. pi8| , (6.10) plays the role of an irrelevance criterion for 
weak, short-range correlated randomness that couples to the surface energy density (and is restricted to the surface). 
If it is satisfied, the fixed point describing the special transition of the pure model is stable with respect to this kind 
of randomness., so that such random "surface-enhancement disorder" should be irrelevant at the special transition. 
According to our numerical estimates, this irrelevance should indeed apply. It has been verified by Monte Carlo 
simulations recently Em. 



VII. ORDINARY TRANSITION 



In our analysis of the asymptotic critical behavior at the special transition it turned out to be advantageous to set 
the bare and renormalized surface enhancements to their respective critical values cq = c'^ and c — 0. The benefit 
was that we did not have to deal with renormalization functions depending on two mass parameters m and c, a fact 
which facilitated the computation of the required Feynman graphs considerably. 

In the case of the ordinary transition we must study the limit c/m — > 00. For the sake of achieving a similar 
simplification, it would be desirable to set c — 00 (or cq = 00) from the outset. In doing so one is faced with a known 
difficulty: Studying the functions G*-^'^^^ with co = 00 does not easily give access to surface critical exponents via 
the RG equations of their renormalized analogs because these bare functions as well as the renormalized ones with 
c = 00 satisfy Dirichlet boundary conditions. Fortunately it is known from previous studies based on alternative RG 
approaches |jlJ,|l5|,|l8|,|o| how this problem can be overcome: one must study the functions 





N 




" M 


C?W^^)(a;i,...,rM)^^ 


n^'^^c^,) 












.fe=i 
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where (9„ means the derivative along the inner normal. The functions ^(^•^^) with M > do not vanish for cq = oo, 
and the scaling dimension of dn4> yields the sole missing surface exponent .p] 

That the relevant information can be obtained from these functions can be seen either by expanding the bare 
cumulants G^^'*'^) in powers of Cq ^ or by noting that because of the Dirichlet boundary condition dnip is the leading 
operator appearing in the boundary operator expansion |l5| , [l^ of (j>. 



A. General considerations and the limit c/m — > oo 



Let us denote the functions Q^^'^^ with 



Co 



^{N,M) 

OO as ypo 



Although we shall not present a complete analysis of 



the c-dependent normalization conditions of Sec. IIIB and of the crossover from special to ordinary surface critical 



behavior here, we will at least verify that this renormalization procedure is consistent with the one based on the 
Qx'^^\ a scheme whose results were briefly described in Ref. and which will be exploited below. 

We start by performing the mass renormalization as described in Appendix Al, and introduce a{p]m,CQ 
analog of ctq (p) , via 



the 



(7^°'^) [p; mo(m), uo, cq] = [k + cq - C7{p; m, cq)] 



(7.2) 



Its expansion to order Uq is given in ( A16 ) of Appendix Al 



Assuming that the renormalized surface enhancement c has an arbitrary value < c < oo , we imagine that the 
surface-enhancement renormalization has been carried out in the way explained in Appendix A 2 . Substituting the 
resulting form ( |aT9| ) of [G(°'2) (p)] into ( |5^ ) yields 



[Zi{u,c/m)Z^{u)] 



- lim — — \a(p: m,c + 5c) — afO; to, c + 5c)] 
P-*o p dp ^ ■' 



(7.3) 



We wish to study what happens to the perturbation expansion in u of the right-hand side of (7.3) in the limit 
c/to oo. To this end, we set m — I and let c — s- oo. Then the free pr opa gator — namely (2.6), with cq and kq 
replaced by c and k, respectively — goes over into the Dirichlet propagator (A3). Further, the perturbative corrections 
caused by the shift Sc to the term inside the square brackets of ( |7.3| ) vanish as c — + oo. Hence we have 



lim \(t{p; m,c + 6c) — a{Q; m,c + dc)] = a^^^ {p; to) 

c — >oo 



.(D) 



(0;to) 



(7.4) 



with a^^'^{p;m) — (T(p;m, oo). As we have seen, the graphs of (7^^-'(p;m) are obtained from those of a{p;m,c) by 
associating with all full lines t he D irichlet propagator Gd rat her t han the c-dependent one ( |2.6| ), and with all full 
lines labeled "s" (cf. Appendix Al) its surface part, given by ( A12| ) with cq = oo. Note that these graphs are not in 
general uv finite at d = 3. But subtraction of their values at p = 0, which is provided by the last term in (7.4), is 
sufficient to make them so. In other words, in the limit c — > oo, surface-enhancement renormalization reduces to an 
additive renormalization. 

To see how this relates to our approach based on the cq = oo functions Gx'^^K we return to the representation 
( ^.5|) of (To in terms of the self-energy S. Since the denominator of the fraction in (|]^) becomes one for cq = oo, we 
have (7^^^ = g^T[Gr,]g, where T[G] is the T-matrix introduced in (5^). Now the reduced propagator ( ^ ) can be 
written as 



d - 

Kp;z')^e ^—Gb{p;z,z') 
oz 



(7.5) 



Thus we get 



g'^^'^\p-m^{m)] = -n + a(''\p-m) 



(7.6) 



dime nsion A[ei] of the surface energy density £i at the ordinary transition is exactly given by A[ei 
: a - 1 and alf = a - 2 



d, 

respectively I6C|,|15|. The other surface exponents are 



^"Since the sc alin ; 
the analogs of (p. 3: 



and (|6.3b[) re ad al"^ 
given by the scaling relations (6.1). 

'^^A simple way to see this is to note that such corrections involve free propagators with points on the surface. Dimensional 
arguments lead to the same conclusion. 
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where it should be remembered | |T7| , p^ , |60t that dzdz>Gj:,ip', z, z') has a contribution of the form [—S{z — z')]; we 
have dropped the impHed singularity [—(5(0)] in the zero-loop term {—n), interpreting dnG^dn as the limit of 
dzdziGi:,{p; z, z') as z, z' ^ with z ^ z' . Combining these findings with (7.3) and (7.4) yields 



[Zi(u,oo)Z,(u)]-i = - Ihn ^1 [^(^■2)(p) -e^)(0) 



(7.7) 



Next, let us recapitulate our renormalization scheme for the q!^'^^^ | ]38[ |. Aside from the previous bulk renormal- 
ization functions, it involves a renormalization factor Zi^oa{u), which enters the definition of the renormalized surface 
operator: 



-{N+M)/2 



(9„0)rcn = [^1,00^0] ^^'^ dnCj) , 

and of the renormalized functions: 

One evident normalization condition is 
The other, 



Qp2 ^oo,rcnV-Pi '"-i "/ 



1 



(7.1 



(7.9) 



(7.10a) 



(7.10b) 



0=0 2ra 

(suggested by the corresponding zero-loop result) serves to fix Zi^ao- In conjunction with (|7.9|) it implies the relation: 



P^O p Op 1 



whose comparison with (7.7) reveals that 



Zi{u,(X))Z^{u) 



lim 

c/m — ^00 



Zi{u,c/m)Z^{u) ^ Zi^oq{u)Z^{u) 



(7.11) 



(7.12) 



to any order of perturbation theory. 

We introduce the analog of the exponent function 77^^ by 



?7i,oo(u) = fjiooiu) = 

om 



lnZi,oo(M) = P{u) 



9 In Zi^ao{u) 

du 



(7.13) 



where u is the rescaled coupling constant of ( |5.8| ). The fixed-point value of this function, ?7i,oo('u*), is related to r]°^^ 
via Hi (cf. Ref. |l|) 

^r' = 2 + ^r(«*)+^0(^*), (7.14) 
as we shall verify below. Reasoning in a standard fashion, we find that 

as rn —> (or u — > u*), with fixed bare interaction constant uq (and cq — 00). 



The renormalized functions Gcx>,'rcn satisfy the analog of the CS equation ( [4.3a| ): 

G[^'tiH;m,u)^AG.^ 



m- h P(u)— H r]^[u} + — ?7i,oo(w) 



dm 



du 



(7.16) 



in which the inhomogeneous term AGco,rcn is defined just as AGion in (1-31:), but with Gr^'.sp replaced by 



G^ren^'°\ the Corresponding cumulant with an insertion of 5 /y </>^. Neglecting AGoo.ron, we can exploit in the usual 
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fashion the resulting homogeneous CS equation together with the asymptotic forms (4.16) and (7.1£) of and Zi.oo 
to conclude that the bare cumulants behave as 



[mx, mr) 



(7.17) 



near criticality. That these scaling forms carry over to the asymptotic behavior of the functions G^^'*^^ {x, r; mo, uo, cq) 
near the ordinary transition can be seen in the ways mentioned in the introduction to this section and expounded 
in Refs. (use either the expansion of the bare functions in powers of I/cq or the boundary operator expansion). 
Here we shall present an alternative derivation, which is based directly on our c-dependent renormalization scheme. 

First, we need the asymptotic scale dependence of the variable c = c( m) n ear the ordinary fixed point. This can be 
conveniently obtained from the reformulated normalization condition ( ^.la ). The bare function G^^'^^{p = 0) = xii 
appro aches a finite value Xii'^(woj cq, A) as T — > Tc [m — > 0) with fixed uo an d cq < Cq . Using the limiting behavior 
( 7.12 ) of Z^Zi for c/m — + oo together with the asymptotic forms ( 4.1(j ) and ( [7.15 ) of and Zi^^o, we arrive at the 
relation 



(c + m) m^+^i'-f"*) - x^i"^ , 



which yields 



The second ingredient we shall need is the asymptotic behavior of the dimensionless function 
^-iN+M)[d-2),2 q[nm) (a.^ ^. c) = GL^/^) [mx, mr; 1, u, c/m) 



(7.18) 
(7.19) 
(7.20) 



as c = c/m oo. Based on our knowledge of the 1/co expansion (cf. the analogous considerations in Sec. Ill C 6 of 
Ref. [0), we anticipate that 



G[^/^'\x,r;l,u,c) 



(7.21) 



where ri2 = ri — r2- When these results are inserted into G*-^'^'^-* = Z^^^'^^^^ zf'^^^Gl:^n'^''\ each one of the M 
surface operators (f)s is found to contribute a factor 



to the prefactor of Hence we recover indeed the familiar scaling form [cf. ( [7.17 )]: 

G^'''"\x,r) ^m(''P+''^°'"^/'' J^(^^''Hmx,mr,u*) . 



(7.22) 



(7.23) 



In the special case {N, M) = ( 0, 2), a contribution m*-''ii 
which we have suppressed in (7.23). 



TZ{u*) S{ri2) and similar subleading ones oc S{ri2) appear. 



B. Results of perturbation theory to two-loop order 



We now turn to the explicit calculation of the renormalization factor Zi,oo(u).^0('u) up to two-loop order, restricting 
ourselves again to the case d = 3. 

Setting Co = oo in (A16) gives us the perturbation expansion of a'-^^ to second order in uo- This we insert into 

( [7.6[ ), and the so-obtained form of G^(x>''^^ then into (7.11). There are two simplifying features we can benefit from. 
First, the one-loop graph of a^^^ differs from its c = counterpart by a minus sign. This means that the term linear 
in uq agrees with its counterpart for Zl^Z^. Second, as we shall show in Appendix the contributions from the 
two-loop graphs (3) and (4) of Fig. 3 cancel. Hence we get 



Zl ooZd, — 1 



2 Mo 



12 87rm 



2m 



_9_ 



p=0 




1 71 + 2 

2^ 18 



il [/2(m2)-m2/3(m2)] . 



(7.24) 



23 



The required Mg term is easily calculated (see Appendix One finds 



with 



107 7 4 

C= In- - 0.094299 = -0.105063 

162 3 3 



(7.25a) 



(7.25b) 



Upon expressing uq in terms of the rescaled renormalized coupling constant u~u (ri + 8)/48m7r [cf. ( [3.§| )], the result 
becomes 



7 ( \7 ( \ , , n + 2 , ^ 12(n + 2) , n + 

^l,oo(w)Z0(u) = 1 + — ^ U + — ^TT^ I C 



2(n + 8) (n + i 
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(7.26) 



From it the exponent function appearing on the right-hand side of (7.14) can be deduced in a straightforward fashion. 
One obtains 



r;,, (u)-,7|| {u)-2 + g) " (n + 8)^ + 96 



(7.27) 



The correspondi ng se ries expansions of the surface e xpo nents A^"''^, ij'^'^, PT'^, 7ii'*: iT'^: ^T'^j a-nd 6 °]'^ f ollow again 
by substituting (7.27) together with the expansions (|6.2|) of and into the scaling-law expressions (6.1). 



C. Numerical estimates for the surface critical exponents of the ordinary transition 



Following the strategy described in Sec. one can analyze the above power series for th e critical exponents of the 
ordinary transition and extract numerical estimates. The results are shown in Tables VIITX], where the entries have 
the same meaning as in Tables (Sec. 0). As before, the fixed-point values u*{n) of Refs. [|l) and obtained 
by Pade-Borel resummation of the two-loop result for the /3 function, were used. 



TABLE VIII. 
u* = 1.632. 



Surface critical exponents of the ordinary transition for d = 3 and n — 0. As fixed-point value we used 





O1/O2 




[0/0] 


[1/0] 


[0/1] 


[2/0] 


[0/2] 


[1/1] 




'711 


2.5 


2.0 


2.00 


1.796 


1.815 


1.715 


1.734 


1.660 


1.660 


Ai 


4.4 


7.8 


0.25 


0.352 


0.364 


0.375 


0.380 


0.382 


0.394 


V-L 


3.6 


2.6 


1.00 


0.898 


0.907 


0.870 


0.877 


0.859 


0.843 


Pi 


-1.9 


-1.6 


0.75 


0.852 


0.864 


0.799 


0.790 


0.817 


0.782 


7ii 


0.0 


0.0 


-0.50 


-0.500 


-0.500 


-0.424 


-0.434 




-0.388 


71 


15.4 


-11.4 


0.50 


0.653 


0.681 


0.663 


0.662 


0.664 


0.680 




2.5 


3.1 


1.67 


1.780 


1.788 


1.825 


1.832 


1.854 


1.870 


Sii 


2.1 


2.7 


0.33 


0.424 


0.433 


0.466 


0.476 


0.504 


0.504 



TABLE IX. Surface critical exponents of the ordinary transition for d = 3 and n — 1. As fixed-point value we used 
u* = 1.597. 





01/02 


Ou/02^ 


[0/0] 


[1/0] 


[0/1] 


[2/0] 


[0/2] 


[1/1] 




'?ll 


2.3 


1.8 


2.00 


1.734 


1.765 


1.618 


1.655 


1.528 


1.528 


Ai 


3.0 


5.0 


0.25 


0.383 


0.404 


0.427 


0.440 


0.450 


0.464 


Vi- 


3.0 


2.2 


1.00 


0.867 


0.883 


0.823 


0.837 


0.801 


0.779 


/3i 


-2.5 


-1.9 


0.75 


0.883 


0.904 


0.829 


0.815 


0.845 


0.796 


711 


0.0 


0.0 


-0.50 


-0.500 


-0.500 


-0.402 


-0.418 




-0.333 


71 


5.7 


-40.4 


0.50 


0.700 


0.749 


0.735 


0.742 


0.742 


0.769 


Si 


2.2 


2.7 


1.67 


1.815 


1.829 


1.883 


1.898 


1.941 


1.966 


Sii 


1.9 


2.5 


0.33 


0.452 


0.468 


0.514 


0.533 


0.582 


0.582 
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TABLE X. Surface critical exponents of the ordinary transition for d — 3 and n — 2. As fixed-point value we used u* = 1.558. 





O1/O2 


Ou/02^ 


[0/0] 


[1/0] 


[0/1] 


[2/0] 


[0/2] 


[1/1] 






2.2 


1.6 


2.00 


1.688 


1.730 


1.545 


1.598 


1.422 


1.422 


Ai 


2.4 


3.9 


0.25 


0.406 


0.435 


0.470 


0.493 


0.514 


0.528 


Vi- 


2.7 


1.9 


1.00 


0.844 


0.865 


0.787 


0.808 


0.753 


0.727 


Pi 


-3.2 


-2.1 


0.75 


0.906 


0.935 


0.856 


0.840 


0.868 


0.810 


7ii 


0.0 


0.0 


-0.50 


-0.500 


-0.500 


-0.387 


-0.408 




-0.282 


71 


3.9 


42.1 


0.50 


0.734 


0.805 


0.794 


0.814 


0.815 


0.851 


6i 


2.0 


2.5 


1.67 


1.840 


1.860 


1.928 


1.952 


2.019 


2.051 


Sii 


1.8 


2.3 


0.33 


0.472 


0.494 


0.550 


0.579 


0.651 


0.652 



TABLE XI. Surface critical exponents of the ordinary transition for d — 3 and n — 3. As fixed-point value we used 
u* = 1.521. 





01/02 


Ol^/02^ 


[0/0] 


[1/0] 


[0/1] 


[2/0] 


[0/2] 


[1/1] 




m 


2.1 


1.5 


2.00 


1.654 


1.705 


1.489 


1.556 


1.338 


1.338 


Ai 


2.1 


3.4 


0.25 


0.423 


0.459 


0.504 


0.538 


0.574 


0.586 


Vi- 


2.5 


1.8 


1.00 


0.827 


0.853 


0.759 


0.787 


0.714 


0.685 


Pi 


-4.1 


-2.4 


0.75 


0.923 


0.959 


0.880 


0.862 


0.889 


0.824 


711 


0.0 


0.0 


-0.50 


-0.500 


-0.500 


-0.377 


-0.401 




-0.238 


71 


3.1 


16.3 


0.50 


0.759 


0.850 


0.842 


0.880 


0.882 


0.927 


Si 


1.8 


2.3 


1.67 


1.859 


1.884 


1.963 


1.995 


2.088 


2.124 


Sii 


1.7 


2.3 


0.33 


0.487 


0.515 


0.578 


0.617 


0.711 


0.711 



For most of the obtained truncated series expansions, the coefScients do not alter in sign, and the truncated series 
of their reciprocal (i.e., their 'inverse series') display a similar behavior. Of this kind are the series for r]°i^'^ , A™'*, 
j^ord^ ^ord^ jord .^-^j^ ^ _ Q, . . . , 3, and for 7"'^ with n — 2 and 3. Let S[p/,j], with p -I- q < 2, be the values 
resulting f rom Fad e approximants of type [plq\ (and listed in the columns marked [p/q\), and let Sp = S[p/o]- Looking 



at Tables VIII - KI one realizes that the sequences of values S[p/q] associated with each one of these critical indices 
have the following feature: The values move away from sq such that the second-order approximants [p/(2 — p)] give 
values farther away from sq than the first-order ones [p/(l — p)] and that furthermore S[i/i] is the most distant one. 
In other words, either they increase according to 

So < si < S2 < S[i/i] and so < S[o/i] < S[o/2] < ^[i/i] (7-28) 

or else they decrease in the corresponding fashion. In most cases even the stronger chain of inequalities 

So < min{si,S[o/i]} < max {si, S[o/i] } < min {ss, S[o/2] } < max {s2, S[o/2] } < S[i/i] (7.29) 

or its decreasing analog applies. 

The value s^im always comes last in these sequences. Using it to extrapolate the series amounts to anticipating 
that the next (thus far unknown) terms of the power series expansion in u have coefficients of the same sign. This 
assumption might well be true for some of the series (an example of this kind is the bulk exponent ?/), and in view of 
the just mentioned feature of the S[p/g] with p -t- q < 2 it seems legitimate to us to accept it. Accordingly we consider 
S[i/i] to be the best among all those estimates S[p/^] with p + q < 2 for a given exponent that we obtained from its 
individual series expansion. 

For 7™'^ with n = and 1, only the first chain of inequalities of ( |7.28| ) holds. Its inverse series has first-order 
and second-order corrections of different signs, and hence may be treated by the Pade-Borel method. The resulting 
resummation values (the analog of the ones denoted R'^ in Tables agree with S[i/i] up to three decimals. 

In the case of (3°^^ (with n = 0,1,2,3) both the direct and the inverse series are alternating. The results of our 
resummations differ from the values of the [1 / 1] approximants only in the third decimal. (Therefore we have not listed 
them separately.) The series for have zero first-order corrections and hence are not well adapted for estimating 
this critical exponent. 



In order to gain further improved estimates, we follow a similar strategy as we did in Sec. VI when estimating 



the critical exponents of the special transition: we try to exploit the above results in conjunction with the available 
high-precision estimates for the bulk exponents v and r}. To this end we substitute our [1/1] values for ?7||'''^, together 
with the estimates taken from Ref. H, = 0.588, -q = 0.027 (for n = 0), = 0.630, 77 = 0.031 (for n = l),u = 0.669, 
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ry = 0.033 (for n = 2), and v = 0.705, ry = 0.033 (fo r n = 3), into the scaling-law expressions ( |6.l| ). The results are 
given as /(?7||, i^, f]) in the last row of Tables VIII -XI . As one sees, in t hose cases in which the Pade values S[p/q] move 
away from sq in a given direction such that either (7.28) — or even (7.29) — or else their corresponding decreasing 
analogs hold, the estimates /(?7||, J^, rj) turn out to be displaced even further in the same direction. 

We consider our estimates /(?7||, i^iV) the best we could attain from the available knowledge on the series expan- 
sions, within the present approximation scheme. In some cases they differ significantly from the zcroth-order values sq 
we started from. Like in the case of the special transition, our best estimates agree reasonably well both with the earlier 
ones based on the e expansi on |l^ , |6l|| as well as with more recent computer-simulation results |39|j43| , |59| , |6^ -^ . The 
latter are gathered in Table XII . For references to earlier numerical estimates and their comparison with e-expansion 
results, the reader may consult Ref. [|l8|. 



TABLE XII. Monte Carlo estimates for the surface critical exponents of the ordinary transition in d = 3 dimensions. 

n = 



711 


-0.38(2) 


Meirovitch & Livne, 1988 




-0.353(17) 


De'Bell et al, 1990 H 




-0.383(5) 


Hegger & Grassberger, 1994 


m 


71 


-K0.718(8) 


De'Bell & Lookman, 1985 E 






-f0.687(5) 


Meirovitch & Livne, 1988 m 






+0.694(4) 


De'Bell et al, 1990 ||6gl 




+0.679(2) 


Hegger & Grassberger, 1994 ji^] 



n = 1 



/3i 


+0.79(2) 


Kikuchi & Okabe, 1985 H 






+0.78(2) 


Landau & Binder, 1990 lo 






+0.75(2) 


Ruge et al, 1993 @ 




+0.807(4) 


Ruge & Wagner, 1995 [El 




+0.80 + 0.01 


Pleimling & Seike, 1998^ 




71 


+0.78(6) 


Landau & Binder, 1990jp 






+0.760(4) 


Ruge & Wagner, 1995 [MT 




+0.78 ± 0.05 


Pleimhng & Selke, 1998^ 


1 


7ii 


-0.25 + 0.1 


PleimUng & Selke, 1998 W 


1 


<5i 


+2.00(8) 


Kikuchi & Okabe, 1985 @ 





n = 2 



/3i 
71 



-0.84 

^ 2 
- 3 



Landau et al., 1989 
Landau et al, 1989 



|35| 
1 351 



Specifically, our estimates j°l'^{n — 0) ~ —0.388 and 7™'^(n = 0) ~ 0.680 for the polymer universality class (n — 0) 
are in excellent agreement with the recent (apparently very precise) Monte Carlo estimates 75'5^'^(n = 0) = —0.383(5) 
and 7°'''^(n = 0) = 0.679(2) by Hegger and Grassberger [Q. Likewise for the Ising universality class, our numerical 
values /3°'-d(n= 1) ~ 0.80 and -f1"^{n = l)~ 0.77 are very close to the Monte Carlo estimates PT'^^ 1) = 0.807(4) and 
7°'-d(n=l) = 0.760(4) of Ruge et al. Q. Landau and Binder's earlier Monte Carlo estimates @ /3°'-d(n=l) ~ 0.78 
and 7°'''^(n=l) = 0.78(6) are slightly smaller and larger, respectively. The more recent ones by Pleimling and Selke 
p9[ coincide within their error bars with those of Ref. [Q and our best estimate. 

There also exist some experimental results^ with which these theoretical Ising values can be compared. Sigl and 
Fenzl ||6^] were able to extract the value /3i = 0.83 ±0.05 from capillary-rise experiments on the transition from partial 
to complete wetting in critical mixtures of lutidine and water with different amounts of dissolved potassium chloride. 
Using the technique of x-ray scattering at grazing incidence |6^j6^,|7^ , Mailander et al. investigated the 
surface critical behavior of a FeAl alloy at its B2-D03 disorder-order transition |72| - [7^ . The values r^n — 1.52 ± 0.04, 
f3i = 0.75 ± 0.06, and 711 = —0.33 ± 0.06 they found are consistent with our estimates^ ?7||''^(n = 1) ~ 1.53, 
/31"^{n=l) ~ 0.80, and '^'lf{n = l) ~ -0.33 (taken from the last column of Table S). 



'^^For background and a review of experimental work on surface critical behavior, see Ref. | |6q |. 

^^The case of the B2-D03 transition is more complicated than that of the A2-B2 transition, for the D03 structure involves 
four sublattices and hence a larger number of composition variables |Q. Nevertheless the B2-D03 transition is expected to 
belong to the Ising universality class | j73|j74| ; see the note added in proof in Ref. [Q. 
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An x-ray scattering experiment has also been performed on the A2-B2 disorder-order transition in a semi-infinite 
FeCo alloy that is bounded by a (001) surface ||7^]. This yielded f3i = 0.79 ± 0.10, in conformity with the above 
theoretical values for (]°'^'^{n = 1). Yet it should be noted that the chosen (001) surface breaks the symmetry of 
interchanging the two sublattices |7^-|7^ . Therefore the Hamiltonian one encounters in a coarse-grained continuum 
description of the large-scale physics is not invariant under a sign change (/> — > — of the order parameter and will 
generically have surface contributions involving odd powers of <j) and its derivatives ||79| , [l8|] . In particular, surface 
contributions linear in 0, i.e., a surface ordering field gi ^ 0, normally should be present, and since gi is a relevant 



scaling field, the asymptotic critical behavior must be characteristic of the norma rather than the ordinary 

transition ||7^,|7||. 

In their experiment, Krimmel et al. |75|| actually found evidence of the presence of such a surface ordering field gi. 
On the other hand, they did not observe the crossover to the normal surface transition. The reason seems to be that 
gi is rather small. In order to see clear manifestations of this crossover or even verify the true asymptotic behavior, 
one must therefore resolve a temperature regime fairly close to Tc. The one studied in the experiment was apparently 
not close enough, a possibility which has already been suggested by the experimentalists themselves A recent 

reanalysis of their data indicates that these are even better consistent with the behavior one should expect near 
Tc when the scaling variable (7i|t|~^i is still small (so that the crossover to the normal surface transition has not 
yet set in) than the original analysis by Krimmel et al. revealed. 

The experiments | |7l| on the B2-D03 transition of FeAl also require a comment. Just as in the measurements on 
FeCo Izsj] , a small amount of long-range order near the surface was found to persist at and above Tc. It is tempting to 
attribute this again to the presence of a surface ordering field gi (cf. Ref. [p6|). However, the orientation of the surface 



plane of the FeAl crystal investigated in Ref. |7l| was symmetry preserving in the sense of Refs. |76 and ||7 



surface contributions breaking the (f> — > —(f> symmetry of the Hamiltonian should not occur. Thus, if the explanation 
of the experimental findings must indeed be sought in the presence of a surface ordering field, then the question of its 
origin arises.]^ It appears that further theoretical and experimental work is required to clarify this issue. 

X-ray scattering experiments have also been performed on a NH4Br single crystal [Q. The authors argue that 
the critical fluctuations at the observed order-disorder transition should be described by the three-dimensional Ising 
model, but also point out that the transition is coupled to a first-order displacive transition. The effective exponents 
they measured, ?7|| — 1.3 ± 0.15 and /3i = 0.8 ± 0.1, are compatible with the theoretical predictions for the n — 1 
ordinary transition. In view of the coupling to the displacive transition it is however not clear to us how serious such 
a comparison can be taken. 

Our estimates for n = 2 and 3, given in Tables ^ and XI, also conform nicely with the previous e-expansion 
estimates gathered in Table VI (p. 186) of Ref. from which we quote the value 77||"^(n = 2) ~ 1.38 as an example 
(to be compared with our present best estimate ~ 1.42). For n = 2, there exist some recent Monte-Carlo results by 



Landau et al. [|65[ , as mentioned in Table XII . For a comparison with series-expansion estimates for the cases n = 2 
and 3, we refer to Table VII of Ref. ||l|| and the original work [831. 

Our values 7yjj"^'^(n = 3) ~ 1.34 and /?5"'^(n = 3) ~ 0.82 are fairly close to the estimates Tj'^^''^in = 3) ~ 1.29 ± 0.02 

and /35"^'^(n = 3) ~ 0.84 ± 0.01 Diehl and Niisser |6l[] obtained from Fade approximants that exploited the results 
of both the e expansion and the d — 2 expansion to second order. We are not aware of any recent Monte-Carlo 
predictions for surface critical exponents of the n — 3 ordinary transition. On the experimental side, there is the 
result (3i = 0.825lJ^;j^4g Alvarado et al. found for a Ni(lOO) surface using spin-polarized low-energy electron 
diffraction. 



VIII. CONCLUDING REMARKS 



In this work we have extended the massive field-theory approach for studying critical behavior in a fixed space 
dimension below the upper critical dimension to systems with surfaces. We have carried out two- loop calculations 
for the ordinary and special surface transitions in d = 3 bulk dimensions and performed a Pade-Borel analysis of the 
resulting series for the respective surface critical exponents. The behavior of the truncated series we have obtained 



^''A justification of the name "normal surface transition" may be found, for example, in Refs. [£0| and | |l9| . 

^^Any real surface will, of course, not be ideally planar. Hence the symmetry invoked in proving the absence of symmetry- 
breaking terms in the Hamiltonian will not be strictly realized. Nevertheless, one would not expect such unavoidable 
imperfections to manifest themselves through symmetry-breaking contributions proportional to the surface area, unless the 
crystal was not carefully prepared and its surface not well-defined. 
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and analyzed, though less good-natured for some than for other exponents, is in general very similar to what one 
finds for those of bulk exponents at the same two-loop order of truncation. We take this as a clear indication of the 
potential power of the approach: whe n pus hed to an order of perturbation theory that is comparable to what has 
been achieved for the bulk exponents [[3y2l|- |2^ , |85[] and investigated by the same sophisticated techniques based on 
Borel summation and large-order analysis, it should yield similarly precise numerical estimates. 

One motivation for the present study was to see whether the field-theory results might be reconciled with the small 
values of ~ 0.5 found in recent Monte Carlo simulations for the crossover exponents pjd = 3, n) with n — 



and n = 1. Our present best estimates $(3,n = 0) ~ 0.52 and $(3,n=l) ~ 0.54 (cf. Tables III and IV) are indeed 
much lower than the original ones based on the e expansion (which were ~ 0.67 and ~ 0.68, respectively [|l6|,|l8|), and 
as we have seen, a Pade-Borel analysis of the e expansion to order yields comparatively low d = 3 estimates. That 
the original e-expansion estimates for $ were ~ 20% greater than our present ones seems to be due to the unusual 
largeness of its 0{e^) terms, which entails that the value of the truncated power series at e = 1 gives a rather poor 
approximation for ^(d — S). This problem exist, of course, also for the other surface exponents that derive from the 



same RG function t]c as $ (such as a'^, cf. Tables HI and IV). For the remaining surface exponents of both transitions, 
the O(e^) terms are much smaller, so the values of the truncated series at d = 3 turn out to be much closer to our 
best estimates. 

In those cases in which there is little difference between the e-expansion values given in Ref. ||l8| and our best 
estimates here, one may say that these field-theory values have been put on a more reliable basis by our present 
analysis because of our use of better extrapolation procedures based on Pade-Borel summation techniques. 

To give error bars for our estimates of surface critical exponent is a rather delicate matter. If we took as a measure 
of uncertainty for the value of any given one of them the spread of values of the various extrapolations of the 0{u'^) 
series expansion, then a reasonable guess might be a typical accuracy of a few, say, 5%. What appears to be needed 
most for an improvement of the accuracy and more reliable error bars is the computation of the series coefficients of 
the surface exponents to a higher order in perturbation theory. 

There is an additional problem one is faced with in massive field-theory approaches to systems with boundaries that 
should be mentioned: the appearance of further mass scales such as the renormalized surface enhancement c. Having 
to deal with more than one mass parameter, namely with m and the ratio c/m, makes calculations rather cumbersome. 
Fortunately, we have found ways to study directly the asymptotic cases c/m = and c/m oo corresponding to 
the special and ordinary transitions, respectively. Hence one gets back to single-mass problems. Nevertheless the 
technical problems that must be overcome to extend the calculations to higher orders of the loop expansion require 
considerably more effort than in the bulk case. 

It is our hope that the present work might serve as a useful basis and starting point for further analyses that 
ultimately could lead to quantitiative field-theory results for surface critical exponents and other universal quantities 
of a precision as good as in the bulk case. Finally, we would also like to express our hope that our work might spur 
further experimental work as well as simulations, the latter especially for higher spin dimensionalities. 
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APPENDIX A: TWO-LOOP CALCULATIONS FOR THE SPECIAL TRANSITION 

In this appendix some details of our two-loop calculations for the special transition will be presented. 

1. Mass renormalization 

Our starting point is the Feynman graph expansion of [G^^'^\p)]^^ to two- loop order, as given by (|5.3|) and (^.7|). 



A useful observation is that the term C3(p) in (5.7) can be combined with the one oc [Ci(p)]^ as 




^ f ^ [Clip)? \ I \ I 

K0+ Co - 
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where the Hne labeled "D" represents the difference 



G{p;z,z') 



Co + Kq 



^-Ko(z + z') 



which is nothing but the Dirichlet propagator 

G'd(p;^,^') = 



2ko 



g-Kolz-z'l _ g-Ko(z + z') 



(A2) 



(A3) 



We now perform the mass renormalization. That is, we express the bare mass mp in terms of m and uq. For our 
choice of bulk normalization conditions made in Sec. Ill A, one has (see, e.g., Ref. pl|) 



where 



mn = m 



n + 2 
6 



uq h (m 



n + 2 
18 



ul [h{m')^m'h{7n')\ +0{ul) 



1 



(A4) 



(A5) 



and 



with 



/2(m2) = /2(0,m2) 



1 



/fei ifc^ {kl + to2)(A:| + m2) [{ki + fca + fc)^ + m^)] 
are familiar bulk Feynman integrals. Unlike Ii and /2, the integral /a is finite in three dimensions. 



1 
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The factor l/87rm may be identified as the value for d = 3 of the one- loop integral 

1 



£)(m) 



(fc2+m2)2 



(A6) 



(A7) 



(AS) 



(A9) 



(AlO) 



of the bulk vertex function f [^^'[^.({0}). In the massive field-theory approach, it is common to absorb this factor in a 
properly define d re normalized coupling constant. We did this when introducing the rescaled renormalized coupling 
constant u via ( |5.8| ). ( Not e that the rescaling factor &„(d) in this equation becomes [6/(n -I- 8)]87r when d = 3.) 

So far (namely, in ( |Al[ ) as well as in the ma in tex t) w e used the convention that the lines in Feynman graphs 
represent the free and reduced propagators ( |2.6| ) and (|5.6| ), respectively. For the purpose of determining the results 
of the mass renormalization, it is preferable to use a different one in which the full and dashed propagator lines 
denote the modified expressions one obtains from the previous ones through the replacement ttIq m?, i.e., of kq by 
K = \/ + rn? ■ 

Let us decompose G into its bulk part 



and its surface part 



Gb(p;z,z')-— e--^'' 



4(p;z,.') = ;^V^e-(^+^') 

2k K + Co 



(All) 



(A12) 



writing 
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G{p\ z, z') = Gh{p; z, z') + Gs{p; z, z') , 



(A13) 



and label lines representing Gb and Gg by "b" and "s" , respectively. Since the term linear in uq in ( A4 ) is given by 
the Feynman diagram 



(AM) 

it is clear that the perturbative corrections originating from this term provide subtractions for the tadpole subgraphs 
such that only their surface part 



s 





remains. Thus the self-energy a introduced in (7.2) becomes 



(A15) 



a{p] m, Co) 



(72(p;m, Co) 



with 



0{ul) 



(A16) 



(72(p;m, Co) 




1 n + 2 



2k 18 



(A17) 



The latter combination, a2, is uv finite for d < 4. However, the contributions of the other three graphs in ( [Al(j ), 
which we shall denote as cti, 173, and (T4, respectively, are not because they still contain an uv divergence that must 
be absorbed through the renormalization of the surface enhancement. 



2. Surface-enhancement renormalization 



Upon substituting (5^°'^^ from (7^) into the normalization condition ( 5.1a ), one is led to an equation for the shift 
5c introduced in (|3.6D: 



6c = ^[Zi{c)Z^] ^ — l| (c + m) + (t(0; m, co = c + Sc) 



(A18) 



This can be solved by iteratio n to determine the dependence of co = c + Jc on c and to, order by order of perturbation 
theory in uq (or u). Using it (7.2) can be rewritten as 



[G(°'2) (p)] ^ ^ K - m + [Zi {c)Z^r'^ (c + to) - [a{p; m,c + Sc) - a{0; m,c + Sc)] 



(A19) 



Graphically the result can be interpreted as follows: To obtain [G'^^'^-' (p)] ^ one can set 6c ^ if each graph of a is 
taken to represent the graph itself subtracted by its value at p = and the corresponding subtractions are also made 
for subgraphs involving lower-contributions to a. This is analogous to the familiar graphical interpretation of mass 

renormalization in the bulk vertex function r^'^\q) (see, e.g., |£^). 

Let us be more specific in the special case c = 0. Upon setting c = in (AI8), we see that the O{uo) term of the 
shift Sc can be written as 



fcW ^ [(ZfZ0)-i](i)m + ai(O;m,co = O) . 



(A20) 



The O{uo) term [{Zl^ Zii,)~-^]^^^ will be calculated in the next subsection. The second contribution to 6c'^^^ is uv 
singular for d = 3. It provides the singular part of the counterterm that cancels the uv divergence of <Ti(p; to, 0). To 
show this, let us calculate this quantity. We have 
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Tio n 



Kd-iJ[-;d 
m 



with Kd = 2i-'^7r-'^/Vr(d/2), where 



2 Jo (Vp^+x/^TT) VyTT 



dy L 



Subtracting from the latter integral its value at P = yields the finite d — 3 expression 

1 



[J(P;d)- J(0;d)]d=3 = - In- 



- ^ -ip2 + C'(p4) . 
4 



These results can be substituted into the Dyson form (7.2) of G^'^'^^ to obtain 



+ 2 liQ m + K 

In 

3 Stt 2to 



(A21) 



(A22) 



(A23) 



(A24) 



At two-loop order, the first-order shift Sc^^^ must be taken into account in ai. It produces an 0{uq) contribution 
of the form 



^1 Co = 0) , 

OCq 



(A25) 



whose uv singular part must remove the divergence of the tadpole subgraph ( |A15| ) of 174. Since the same subgraph 
appears (twice) in (T3 one may be surprised at the apparent absence of the corresponding subtraction terms one obtains 
through replacement of one or both of these subgraphs by their singular parts^, and one may wonder whether 
be uv finite. In fact, the analog of 0-3 that has the Dirichlet line "D" replaced by a full propagator line would need 
such subtractions, and those would occur, had we not reorganized the expansion in the manner of (Al) to obtain 
graphs whose internal vertices z and z' are connected by the Dirichlet propagator G'd(p; z, z'). The contribution (T3 



is uv finite bep-use Gd vanishes whenever z or z' approach zero. Replacing one or both tadpoles (A15) in (T3 by their 
singular partsEj oc 5{z) and cx 5{z') would produce vanishing results. 

According to ( |3.2a ), we just need the expansion coefficients of order of [G^'^'^-' (p)]~^ to determine Z^. There- 
fore we shall refrain from deriving the full two-loop resul t for [G*-"'^-* (p)]~^ here. Before turning to the calculation 
of Z^Zij,, let us compute the derivative dai/dco in (A25). This is finite for d — 3, and its 0{p'^) coefficient will be 
needed below. A straightforward calculation gives 



-— (p;m,co = 0) = uq — — - 
oco 12 



uq n + 2 



p' 



■In 



TO 



Mo n + 2 



In 2 



1 

-(l-ln4)i^ 



(A26) 



(A27) 



3. Calculation of Z^^Z^ 



Using the representation (5.2a) of this rcnormalization factor, one can easily deduce from (A21) and ( A23| ) its ©(uo) 
term: 



^^In the dimensionally regularized theory, the singular part of ( |A15| ) has the form (d — 3) ^ &(z)\ if one rather employs a cutoff 
regularization scheme, restricting the parallel momentum integrations by p < A, then the pole gets replaced by In A. 
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dai{p;m,0) 



n + 2 Ufj 



p=0 



12 87rm 



(A28) 



The terms of order Uq can be written as 



1-1(2) /n + 2 Uq \ 1 — ln4 



3 Snm 



2to^ 



i=2 



9p2 



(Tj (p;m, 0) +(5j4 (Ti(0;m, 0) -— (p;to,0) 



(A29) 



p=0 



Th e firs t term is the part of ( |A25D produced by the contribution [(ZiP^'^)-!]^!) of the shift ( |A20| ); the remaining part 
of (A25) is the last term in the square brackets. 
To compute the contribution of (J2 , we start from 



co=0 



n + 2 
18 



00 /'OO 



dz / dz' d'^-^r[GN{x,x')]e'P ' 



where x = {xn,z), x' — {x\.,z'), r = xn — ajn, and Gn denotes the Neumann propagator. Introducing 



we can write the latter as 



i?± = ^r^ + {z±z'y 



GN(a^,a^') = Gb(i?-) + Gb(i?+) 



(A30) 



(A31) 



(A32) 



Since the term ( |A30| ) is not uv finite for d = 3, one must use the form of the bulk propagator for general values of rf, 



Gb(x) = (27r) 
rather than the simpler one for d = 3, 

Gb(x) 



(A33) 



(A34) 



The derivative d/dp^\o may then be performed under the integral sign. Recalling the definition (A17) of &2, one 
thus arrives at 



2m 



9(72 (p; m, 0) 



Qp2 



p=0 



n + 2 
18 



A, /J., p— ± 



where we have introduced the integrals 

rOO roc 

I^^^P^I dz I dz'iz + z')e-"''^'+^'^ I d''-'rGbiRx)G^iR^)G^{Rp) 



and 



00 roo 



J^^''= dz dz'e-^^'+''^ d''-'rr^G^{Rx)Gy,{R^)Gy,{Rp) 



(A35) 



(A36) 



(A37) 



Owing to the bulk singularity contained in the distribution [Gh{R-)]^ , I is not regular at d = 3. Its singularity 
must be canceled by that of the subtrahend oc /2(m^), a term which can be cast in the form 



hjm^ ) _ 1 
2m2 



1 f f 

— dz d'^-\ [Gb(i?)]' 
^ Jo J 



(A38) 



with R — Vr^ + z^. The symmetric double integral dz dz' of / reduces to a single one, if we rewrite it as 
J^^dz^ J^dz^ with z± = z + z' and integrate over This yields 
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him^) 1 



^ dz [(1 + mz)e~"'' -1] J d'^-\ [Gb(i?)]^ 



We may now set d = 3 and perform the r integration. The result is 

, /2(m2) 1 

1 — = z 

2m? 32TT^jn'^ 



with 



dz(e-"^ + mze"™^ - 1) 



-Smz 



3mEi (3mz) 



(A39) 



(A40) 



(A41a) 



where Ei(a:) is the exponential integral function Q. This can be evaluated with the help of the tables and 
One gets 



2 - 71n- 



(A41b) 



The remaining integrals in (A35) are uv finite for d = 3, so one can set d = 3 from the outset. A straightforward 
calculation gives 



J 



5 4 
- + 3 n- 

6 3 



4 3 



and 



31 4 

h 9 In - 

12 3 



(A41c) 
(A41d) 

(A41e) 



The integrals I^^^ and JT±± require a little more effort. Consider / for example. In the coordinates r, 
X = {z — z')/r^ and y — {z + z')/r, the r integration becomes trivial, giving 



1-+ oo 2 2 T — 1-+ 2 

I — S2n m I = m 



OO pOQ 



"'0 



dr dx dy 



OO PX 

dx I dy 

"'0 



with Kx = V 1 + x^. In a similar manner one derives 



OO pX 



dx I dy 



dx 
Jo 



and 



1+ — — 



("^ 


2kx " 






x 




' 2 f 


f ^ 


- 2Ky)^ 




1 




/ 2 ( 

'^x'^y ' 








1 




V 2 ( 




f 2KyY 



(A41f) 



(A41g) 



(A41h) 



(A41i) 



These integrals could be evaluated numerically in the present form. However, by making the familiar Euler sub- 
stitutions, one can reduce the integrands to rational functions. Then the y integrations become standard and can be 
done exactly, so that one is left with single integrations over x. Finally, the latter can be transformed into integrations 
over the bounded interval [0, 1], which is more convenient for numerical calculations. In Table XIII we have listed the 
numerical values of the integrals i^^^ and j^^^ together with those of the analytically known integrals i , j , 
and 
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TABLE XIII. Numerical values of the integrals i^^'' and j^^'' for d = 3. 







H 


- + + 


+ + + 




-0.0137745 


0.1324282 


0.0420126 


0.0239076 




0.0297129 


0.0114128 


0.0075298 


0.0058053 



Substituting the above results and the value (A9) of into (A35) finally yields 



2m 



9(T2(p;m,0) 



71 + 2 / Uq \ 2 



p=0 



3 V87rm 



A 



(A42a) 



with 



A, /I, ± 

17 14 . 

In - + 2^ +1 

162 3 3 



+— -L + j+— + ~ 0.202428 



Next, consider (73. Noting that the integrand is a symmetric function of z and z', we find 



(A42b) 



(A43) 



The integrals one obtains upon performing the derivative d/dp^\^ can be worked out utilizing the mathematical tables 
^ |§. This givesQ 



2m 



p=0 



1 f n + 2 Up y 
4 V 3 Sttto/ 



1\ TT^ 1, 2 

h - In - 

12 2 3 



(A44) 



where Li2(x) is the polylogarithm of order 2 (see, e.g., Ref. [^) with the value Li2(— |) = —0.4484142 . . .. 

It remains to compute the contribution of the combination (74 + ai dui/dcp. The second term implies that from 
the term associated with the lower loop of 172, 



2 ./o 27r 



(A45) 



its value at z' = 0, multiplied by e ^ (the factor produced by the broken lines of iti), gets subtracted. This entails 
that the contribution we are concerned with takes the form 



2m 



_d_ 



p=0 



(T4 + CTi 



dai 
dco 



{p; m, 0) = 



1 / n + 2 Up Y 
4 V 3 87rm J 



94: 



with 
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00 />00 



dz' 



F{z,z') - F{z,0)e 



-2z' 



'0 JO ^ 

where F denotes a sum of exponential integral functions : 

F(z, z') = Ei(2|z - z'\) + Ei(2z + 2z') + 2Ei(|z - z'| + z + z') 
The integrals contributing to g4 can be evaluated analytically. One obtains 



(A46) 



(A47) 



(A48) 



Note that formula 3.1.3.2 on p. 567 of Ref. [p7| is incorrect; its corrected form is given as 3.2.1.2 on p. 617 of Ref. 
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54 = 

Combining the above results yields 

n + 2 Uq 



, 1 3 /I 



1 



n + 2 



12 87rm 



n + 2 (\ 



12 



In 2 + In^ 2 + A 



from which (5.17a) follows upon substitution of 

Mo 6 



87rm n + 8 



(A49) 



(A50) 



(A51) 



4. Calculation Z!? 



We start from the representation (5.2b) of this renormalization factor, perform the mass renormalization, and use 
([7.2[) to express \G^^'^'>\'^ in terms of ct, obtaining 



1 - 



9(t(0; m, cq) 
9co 



(A52) 



c=0 



Upon carrying out the surface-enhancement renormalization, we arrive at 

2 



Z,^ = (2^1 ^0) H In 2 



(7i(0; m, 0) 



3 87rTO 
9^(Ti(0;m, Co) 



3 87rTO 



-ln2 + ln^2 
4 



i=2 



^^^(O; m, Co) 



cfco 



+ 0(Mg) 



(A53) 



co=0 



The explicitly given first C(u§) contribution is the sum of the three terms (9(7i/9co)^, — (Zl^ Z^Y^'' d&i/ dca, and 
— (Z^^Z0)(^'m S^CTi/Sco, whose last one is one part of the contribution proportional to the O{uo) shift Sc'^^^ (the other 
part being the contrib ution propo rtion al to the uv singular tadpole graph ai). The values of (Z^^Z^)*^^^ and d&i/dco 
may be read off from (A5C) and (A_26), respectively; the required second derivative is 



del 



(0; m, co = 0) = — uo 



n + 2 
6 



n + 2 Uq 



p, (k' + m){K') 



3 A-Km? 



(l-ln2) 



(A54) 



The remaining contributions in (A53) can be computed by proceeding along lines similar to those taken in our 
calculation of Z'^^Z^. The contribution of (T2 is given by the derivative d/dco of the graph (A30) at co = 0. It can be 
written as 



9<J2 „^ n + 2 / Mo \ 2 

— (0;m,co = 0) = -^(^j B 



with 



B 



'0 Jo 



6 



oco 



(A55a) 



(A55b) 



where again r = X|| — a;'|. The required derivative {dG^ /dco)co=o — 3(7^ {dG/dco)co=o can be computed in the pz 
representation. Upon setting d = 3 and performing the angular integration, we find 



-^Gix,x') 

OCo 



3-k(z+z') 



' dpp 
2tik? 



Jo{pr) I 



-K,{z-\-z') 



(A56) 



We insert this, together with the decomposition (A32) of Gn, into (A56) and perform the angular part of the 
integration / d^r. This yields 
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B = 2 



with 



/■OO 1*00 7 /"OC /"OO 

Jo Jo "'0 "'0 



-m(_RA+flf,) 



R\R 



(A57) 



(A58) 



where i?A=^ is defined by ( A31 ). 

Changing again to the variables z± — z ± z', we can rewrite the symmetric double integrals Jg dz dz' in k 
and k'^'^ as dz_ dzj^ and dz^ J^^ dz_ , respectively. The inner one of these z integrations can now be 
performed to get 



k^^ = 



00 poo 

dr r 
Jo 



dpp 



Jo{pr) I dz 



-{K+m)z — 2mB. 
^2 



{k + m) ^ for A = — , 
z for A = + , 



(A59) 



where R = Vr^ + z^, as before. Next we make the variable transformation r — > r/z. Then the integrals over z simplify 
to 



1 for A = - , 
z , for A = + , 



standard integrals, which can be found in tables pq ] . In this manner one obtains 

dpp 



00 J poo 



dr r 



«r Jo nl {Kp + l)y/{l + Kp + 2HrY + {PTf 



and 



°° drr f°° dpp 



'^r Jo 



1 ~|~ 2/^-^ 



(1 + Kp + 2K.rY + {prf 



3/2 ' 



(A60) 



(A61a) 



(A61b) 



where as before Kp = yjp^ + 1. To compute these integrals, we had to resort once more to numerical means. Our 
results read 



k— = 0.1538951, 
k++ = 0.0469337. 



(A62) 
(A63) 



The "mixed" term fc + is more complicated since its double z integration does not reduce to a single one. We first 
rescale z, z', and r, so that z/r — > z, z'/r — > z', and mr — > r. This leads us to 



00 poo 



k-+=l drr I ^J^ipr) I dzl dz' e"'' 



Jo 



00 poo 



Jo 



Vk-^'P + i ^{z + z'y + i 



(A64) 



The integral over r is of the form ( A60 ) with A = + and hence can be calculated. Rewriting dz /p dz' as 
Jp°° dy dx with ?/ = z + z' and x = z — z' finally gives the representation we used to determine k ^ by numerical 
integration: 



°° dpp dy dx Ky + Kx + {Kp + l)y 



Jo 



° '^'■^ {[Ky + Kx + {Kp + l)yf + p^^ 



3/2 ■ 



Our result for fc + is 



~ 0.0691008 . 



In conjunction with (A62) it yields 



B = 2 {k— + 2k-+ + k++) ~ 0.678061 . 



(A65) 

(A66) 
(A67) 
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which is the numerical value quoted in (5.2C) 



The other two-loop contributions to can be worked out analytically. Performing first the parallel-momentum 



loop integrations and thereafter the remaining ones over z and z' , one finds thati 



and 



( ^ 

\dco 



f^(0;™,0) 
dco 



^i^)(0;-,o) 



_+2 Up y 

3 8nm I 



U.I -I 



n + 2 



3 Sttto / 



24 



h -ln2 ln3 

12 2 4 



5 , „ 3 , „ 
- In 2 + - In 3 
2 4 



(A68) 



2 MQ y 



- In - + - In^ 2 - Li, 



TT 

12 



(A69) 



In the first version of (A_69), the term in the first pair of squar e bra ckets represents the part of da^/dcp that results 
from the action of the derivative d/dcp on the tadpole subgraph ( |A15|) of (T4 (by which an insertion of (/)g/2 is produced 
there). Whereas this part is uv finite at d = 3, the other one (with an insertion of (j)^/2 at the lower loop) must be 
combined with the remaining term cx cti to obtain the finite result pertaining to the second pair of square brackets. 

Combining these results yields 



Svrm V 4 



3 

n + 2 



VSTrm/ 



A-B 



— flln22-ln2 
3 V4 



(A70) 



from which ( 5.17b[ ) follows by expressing up in terms of u using ( |A5l| ). Just as in the case of Z^^Z^, the number 
Li2(— 1/2) has canceled out. 



APPENDIX B: TWO-LOOP FEYNMAN DIAGRAMS FOR THE ORDINARY TRANSITION 



We first show that the last two two- loop contributions in (AI6), (73(p; m, cg) and <T4{p;m, cq), sum to zero when 
Co — 00: 



a^ip; m, cx)) -I- aiip; m, 00) 
A straightforward calculation gives 



with 



aj (p; m, 00) — — — 



Mo f n + 2 



128 \ 3 / Jp^ Jp,^ K1K2 



■ Lj{k,ki,K2) , j = 3,4, 



nOO rOO 

i3(K,Ki,K2) = - / dz ^-<^+z')'2(n,z+n^z') 

Jo Jo 



4 

K 

1 1 



1 



1 



K K + Kl + K2 \K + Kl K + K2 / K K + Kl K + K2 



11 1 



(Bl) 



(B2) 



(B3) 



Again, the tables of integrals were used; note that formula 2.5.11.11 is misprinted: the argument of the exponential 
function in the integrand should read (c — b)x instead of (6 — cjx. 
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and 



L4(k, Ki,K2) = -- 



«1 Jo 

1 

Adding these contributions yields 

4 



-1 



dz 



' g-2(K2 + K2z') ^g-Kilz-z'l g-Kl(2 + z')^ ' 



1 



1 



K + Kl K + K2 



Kl, K2) 



Kl - K2 



J=3 



(k + + K2)(k + Ki + K2)(k1 + K2) 



(B4) 



(B5) 



Hence the integral /^^ /^^ in (B2) is identically zero because of the antisymmetric form of its integrand. 

There fore the 0(wq) ter m of Z i ,^Z^ is indeed entirely given by the contribution of the quantity (J2(p; m, 00) defined 
in ( |A17| ), as we asserted in ( 7.24 ). To compute it, we can proceed as in our calculation of its cq = analog in Appendix 
^. This leads us to 



2m 



9(72 (p; m, 00) 



p=0 



87rm 



C 



(B6) 



with 



c :^ 3 (r 



h i + r 

162 



(B7) 



Here i are the integrals introduced in (|A41a)-( A41i ). Using the analytic al val ues given in these 

equations together with the numerical ones listed in Table XHI, one arrives at the form (7.25) of our results for 
■Z'l, 00-^0 and C. 
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